{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import shapefile\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import os"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "plt.rcParams['font.sans-serif'] = ['SimHei']\n",
    "plt.rcParams['axes.unicode_minus'] = False"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "def canopy_crow_create(infile,outfile,plotsize = 20):    \n",
    "    fig = plt.figure(figsize=(6,6))\n",
    "    csv_file = os.path.basename(infile)\n",
    "    plotID = csv_file.split('.')[0]  # get from csv\n",
    "    ds = pd.read_csv(infile,encoding='gb2312')\n",
    "    npoints = len(ds)  # get from the csv file\n",
    "    \n",
    "    coods = ds.iloc[:,5:7].values\n",
    "    print(coods)\n",
    "    plt.scatter(coods[:,0],coods[:,1])\n",
    "    theta = np.arange(-np.pi,np.pi,np.pi/50) #必须构建顺时针方向才可以处理\n",
    "    sf = shapefile.Writer(outfile)\n",
    "    sf.autoBalance = 1\n",
    "    sf.field('PlotID','C','10')\n",
    "    sf.field('TreeID','C','10')\n",
    "    sf.field('East_West','N',decimal=8)\n",
    "    sf.field('North_South','N',decimal=8)\n",
    "    \n",
    "    \n",
    "    for i in range(npoints):\n",
    "        # read a,b\n",
    "        a = ds.iloc[i,3]#/2.0 # east-west\n",
    "        b = ds.iloc[i,4]#/2.0  #get from csv, north-south\n",
    "        if a > 10 or b >10:\n",
    "            print(plotID)\n",
    "        treeID = ds.iloc[i,0] # get from csv\n",
    "        crown_x = a*np.cos(-theta)+coods[i,0] #必须构建顺时针方向才可以处理\n",
    "        crown_y = b*np.sin(-theta)+coods[i,1]\n",
    "        plt.plot(crown_x,crown_y)\n",
    "        lcrown_xy = np.column_stack((crown_x,crown_y)).tolist()\n",
    "        sf.poly([lcrown_xy])\n",
    "        sf.record(plotID,treeID,a,b)\n",
    "    plt.xlabel('E(x)')\n",
    "    plt.ylabel('E(y)')\n",
    "    plt.text(-3,23.5,'米')\n",
    "    plt.text(23,-1.5,'米')\n",
    "    plt.savefig(r'./输出数据/树冠投影图.png',dpi=300)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[11.1052999  1.33337  ]\n",
      " [15.5511999  8.4483099]\n",
      " [ 3.22404   12.5551996]\n",
      " [ 9.4088001  1.724    ]\n",
      " [ 1.12738    2.61847  ]\n",
      " [ 9.1027403  4.4907598]\n",
      " [ 7.9523802  7.4126902]\n",
      " [ 1.28083   18.6287994]\n",
      " [16.093399  10.8924999]\n",
      " [ 2.7725999  8.9899397]\n",
      " [12.4624996 17.8083992]\n",
      " [17.6625996  0.558235 ]\n",
      " [ 5.2719898 12.8069   ]\n",
      " [ 2.1186199 17.9269009]\n",
      " [ 7.1174698 11.3796997]\n",
      " [ 5.3183498  2.95105  ]\n",
      " [19.9682999  6.92276  ]\n",
      " [12.1112003 12.1977997]\n",
      " [18.5165997 16.1604996]\n",
      " [17.5067005  1.378    ]\n",
      " [10.7227001 10.9345999]\n",
      " [12.9055004 12.7118998]\n",
      " [11.6405001  0.391265 ]\n",
      " [ 1.2230901  0.53442  ]\n",
      " [13.0902004  6.5057502]\n",
      " [ 7.51195    6.0185699]\n",
      " [15.3094997  6.0210099]\n",
      " [ 6.6721401 14.7740002]\n",
      " [ 0.319936   8.2140102]]\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAGACAYAAACz5SV2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAACI6ElEQVR4nOzddXhcZfbA8e8dn7i7p5F6U3faIi0FSimuC7u47bIsuqzBsvCDZXdZZBd3K9AWKVajSt1SS9rG3TPJuNzfH2lZSuOZyUyS9/M8fZomV86kyZx77/u+50iyLCMIgiAIHVF4OwBBEATBt4lEIQiCIHRKJApBEAShUyJRCIIgCJ0SiUIQBEHolEgUgiAIQqd8LlFIkhQmSdLMEx/fLEnSAyc+XihJ0gjvRicIgjD0+FyiOOFVSZKGAQ2AU5KkFOBZoNmrUQmCIAxBki8uuJMkaTZgBFKBFGAdoJdleZM34xIEQRiKfCpRSJJ0C3A14DrxqUhAB5Se+LcKeFeW5f96ITxBEIQhydcShRZwyLLslCTpHmAhEAK8L8vyP098XZJl2eLNOAVBEIYSlbcD+ClZlq0AkiQlAzcAz9B2VzFXkqSjsix/6c34BEEQhiKfG8yWJCkHWAn8jrZxCoDrgT9IkvS2JEkjvRWbIAjCUORTieJEEvgEuFOW5e9Ofl6W5QZgNnCMtjsMQRAEoZ/42qOng5IkZcuybD/xKQ0nYjzxWOpRrwUnCIIwRPnUYLYgCILge3zq0ZMgCILge0SiEARBEDolEoUgCILQKZEoBEEQhE6JRCEIgiB0yqemx/5URESEnJKS4u0wBEEQBoxdu3bVybLs9rVmPpsoUlJS2Llzp7fDEARBGDAkSSr2xHHFoydBEAShUyJRCIIgCJ0SiUIQBEHolEgUgiAIQqdEohAEQRA6JRKFIAiC0CmRKARBEIROiUQhCIIgdEokCkEQBKFTIlEIgiAInRKJQhAEQeiUSBSCIAhCp3y2KOBQIDscOGpqcNTUYK+txWUw4LJYkC1WZJcThVaLpNWh8NOjiohAFRWFKjoaZUCAt0MXBGEIEYmin7isViwHDmDatRvrkcNYjx7DWlQEdnuPj6WKikI7LB1tRgb6cePQj5+AOjrK/UELgiAgEoVH2YqLaVm3jta16zDv2YN8IimoExLQDhtGwBmzUScltd0pREaiCglB0uuRNFokpQLZasVlteEytuKorcNRW4u9sgLbseNYjx+n8aOlNLz1dtsxExMJmD2bgHlz8Z80CUmj8eZLFwRhEJFkWfZ2DO2aOHGiPBD7UTgaGmj+/HOal6/AmpcHgDYzE/8ZM/CbOAF9Tg6qsDC3nEu227EcOYJp1y5M27Zj/OEHZIsFRUAAQecuIPiiJehzxiFJklvOJwiCb5MkaZcsyxPdflyRKNzDvH8/9W+8Qcuq1eBwoBszhuDzzyNg3jw0CQn9EoPLYsH4ww+0fPsdhu++QzaZ0KSlEXbdtQQvXoxCp+uXOARB8A6RKHyQLMsYN2+h7r//wbxzF4rAQEIuuYSQJRehzcjwamzOViMt335L44cfYsnNRRkaSug1VxP2i1+IwXBBGKREovAx5txcap75B6atW1HFxhL2i+sIueRSlAH+3g7tFLIsY965k/rXXqf1++9RhoYScesthFx5JQoxjiEIg4pIFD7C0dhIzVNP07x8edub7m23EnLFFQPiTdece4CafzyD6YetqJOSiPnjHwmYOcPbYQmC4CaeShRiwV03ybJM07LlFJy7kOYvviD8phtJX/UdYdddNyCSBIB+9CiS33iDxFdfRZIkSm+8kfJ7f4ejvt7boQmC4MNEougGR2Mj5XffTeXDD6NJSyN12adE3XvvgH3WHzBzBqmff0bEnXfS8t13FCy6kNYNG7wdliAIPkokii4Yt2+n8MLFtHy/nqj77yf53XfQZWZ6O6w+U2i1RN55BymffoIqLIzSm2+h6vG/Idts3g5NEAQfIxJFB2RZpuG99yi54Zco/PxI/ehDwn95A5JicH3LdJmZpHzyMaHXXkvjO+9Q8stf4Who8HZYgiD4kMH1rucmssNB1Z/+TPVjfyVg1ixSPvkY3YgR3g7LYxRaLTG/f5i4p5/CnJtL4SWXYMnL93ZYgiD4CJEofsZlsVB2969pWrqU8JtvJuGF5wfsWERPBV9wAcnvvgsOJ8XXXotp9x5vhyQIgg8QieInnK2tlN50M63r1hH9xz8Q9dt7kJRKb4fVr/SjR5HywfuoQkMp+eUvad24ydshCYLgZSJRnOAymym99VZMe/YQ9/TThF11lbdD8hp1fDzJ772LJjWVsjvuwPjDD94OSRAELxKJAnDZbJTdeRfm3XuIf+r/CD7/PG+H5HWqiAiS33gdTUoKpbffgWn3bm+HJAiClwz5RCHLMpUPPohx82ZiH3uUoIULvR2Sz1CGhJD0+muoo6IoveVWrMePezskQRC8YMgnirrnX8Dw1ddE/va3hFx8sbfD8TmqiAiSXn8NSaOh9LbbcTQ2ejskQRD62ZBOFIavv6buhRcIXryY8Jtu9HY4PksdH0/C88/hqKqi/O5f/9iASRCEoWHIJgpbURGVv38EfU4OMY/+RTT36YJfTg6xf30M044d1D73vLfDEQShHw3JROGy2Sj/7b1IajXx/3hmwBT187bgRYsIufQS6l95hdbNm70djiAI/WRIJoraZ5/FcugQsU/8DXVsrLfDGVCiH34Y7bB0Kh54EGdTk7fDEQShH3gkUUiSFCxJ0teSJH0nSdJySZI0kiS9JknSD5IkPeKJc3aXOfcADW+8ScillxI4b543QxmQFHo9cU8/jbOpier/e8rb4QiC0A88dUdxNfAPWZbPAaqAKwClLMvTgDRJkrzSJ1S226l85BFUERFE3X+fN0IYFHTZ2YTf+Cualy+ndZN4BCUIg51HEoUsyy/KsrzqxD8jgWuApSf+/R0w0xPn7Urj++9jzcsj5k9/RBkY6I0QBo2I225Dk5pK9WOPidLkgjDIeXSMQpKkaUAoUAqUn/h0AxDdwfY3S5K0U5KknbW1tW6NxdHYSO0LL+I/cyaBZ57p1mMPRQqtluiHHsRWXEzjBx94OxxBEDzIY4lCkqQw4Dngl0AroD/xpYCOzivL8suyLE+UZXliZGSkW+Ope/E/uFpbxSMnN/KfNQv/GTOofeFFsRBPEAYxTw1ma4CPgYdkWS4GdvG/x01jgSJPnLcj9qoqmj78kJCLlwyK7nS+QpIkou6/H5fBQMPbb3s7HEEQPMRTdxS/AsYDv5ck6XtAAq6VJOkfwGXASg+dt131r7+O7HIRfsut/XnaIUGXlUngOefQ+M67OA0Gb4cjCIIHeGow+z+yLIfKsjznxJ+3gDnAVmCuLMvNnjhvexz19TQt/ZjgRYvQJMT312mHlIhbb8HV2krje+95OxRBEDyg3xbcybLcKMvyUlmWq/rrnABNH3+CbLGIWk4epBsxAv9Zs2h8/wNRB0oQBiGVtwPwJNnppHHpR/hNm4o2La3fz9/U1ER5eTkNDQ2YTCZkWUar1RIcHExMTAzR0dEoB0kHvdArr6Ts9ttpWbOWoAXzvR2OIAhuNKgTRev6DTgqKol+8MF+O6fJZGL37t3s27ePn07xVavVSJKE7SdrDnQ6HSNGjGDSpEnEDvBSIgFnzEYVF0vjhx+KRCEIg8ygThTNn3+OMjycwLlzPX4up9PJ1q1bWb9+PTabjaSkJObPn09SUhIRERFotdoft2tubqa8vJyjR4+Sm5vL7t27GT58OAsWLCA4ONjjsXqCpFQScvHF1D33PPaqKtQxMd4OSRAENxm0icLZaqT1++8JWbIESa326LlaWlpYunQppaWlZGZmMm/ePGI6eKNUKpWEhYURFhbG6NGjMZvNbNu2jU2bNlFYWMiSJUvIHKBTeIMWLqTuuecxfP0N4Tdc7+1wBEFwk0FbPbZ13Tpki4Wghed69DzNzc28/vrrVFVVcfHFF3PVVVd1mCTao9frmTNnDrfddhshISF88MEH5ObmejBiz9GmpqIbMQLDV195OxRBENxoECeKtSgjItCPH++xc9jtdt5//32MRiO/+MUvGD16dK+PFR4ezg033EBSUhLLly+npKTEjZH2n8D587Hk5mKvqfF2KIIguMmgTBSy04lx8xYCZsxAUnjuJa5bt47q6mouueQSEhIS+nw8rVbLFVdcQXBwMMuWLcM+AKeaBsxqW4Bv3LLFy5EIguAugzJRWA4dwtncjP9MzxWpbW5uZuvWreTk5Lh1TEGv17No0SKamprYsWOH247bX7TZ2SjDwzFuFolCEAaLQTmYbdq+HQD/6dM8do7du3fjcrmYPXu224+dmppKcnIyO3bsYNq0aQOqn7ekUOA/dSqmbdu8HcqgYLE72VpQzw8F9RyqMFBcb6LBaMPmcKFWSoQFaEgO82dkXBBT08KZlh6OTj041uYIvmNQJgrzvv2oExNRhYd77Bz5+fkkJycTGhra5bYr9pTz9Ld5VDSZiQvRc9/8LBbndF5OZPTo0Xz55ZfU1dXh7kq6nqYfOxbDypXYq6tRR7dbUV7oQkFtK69vLuSzPRW0WB2olRJZMYGMTQwhIkCDVqXE6nBS32qjsM7IG5uLeGlDAQFaFReMjeNXM1MZFhXg7ZchDBKDM1Hs34+fBwexXS4X1dXVTJvW9R3Lij3lPLQsF7PdCUB5k5mHlrXNauosWSQmJgJQVVU18BLFmLZBffO+fajPOcfL0QwsjUYbT317hI92lKJSKjh/dCwXjItjamo4ek3Hdwon7zy+2FfJst1lfLSjhIvHJ/DgudmEB2j78RUIg9GgSxSOhgYcVVXo+jADqSutra24XC5cLleX2z79bd6PSeIks93J09/mdZooAk904Gttbe1bsF6gHT4cVCosBw8RJBJFt205VsfdH+6lyWTj+ump3DYnncjA7r3J69RK5mRFMScriocXZvPf9cd5Y3MR6/Jq+Mdl45idObAuNgTfMugGs20FBQBoh6V77BxOZ9sbf2M3mvVUNJl79Pmfn2Mg1oJSaLVoEhKwFRZ6O5QB46MdJVzz2jZC/NR8fudM/njBiG4niZ8LD9Dy+/NG8OXdMwn313L9G9t5d2uxmyMWhpLBlyiKigDQpKR47BwhISFotVr8/f273DYuRN+jz5/U0NAAMGBLemhSU3/8vxA69/HOUh74NJdZGZF8dscMRsQFueW42TFBLL9jOnOzonhkxQE+2D4w1+YI3jf4EkVxMajVqOPiPHYOSZJITEykoKAAWZY73fa++VnofzYLRa9Wct/8rE73KzhxZ+SO9RneoElJwVZc3OX3Z6jbVdzAQ8tymZURwcvXTcBf696nwX4aFf+5ZgJzsiJ5ZMUBthXUu/X4wtAw6BKFo6YGVWQEkocf2YwcOZLGxsYf39A7sjgnnieWjCY+RI8ExIfoeWLJ6E7HJxwOB3v27CEtLa1bdy2+SBUdhWy14hJd7zpktjm556N9xIboeP6q8WhVnvmZ1agUPHdlDomhen67dB9Gq8Mj5xEGr8GXKOrqUYVHePw8o0aNIigoiFWrVv04ntCRxTnxbH5wHoVPnsfmB+d1OTV269atGAwGpk+f7s6Q+9XJ/wNHvbiC7chrmwooaTDx1MVjCdZ7tnBloE7N3y8dS3mTmZc3dH5xIwg/N/gSRX29R9dPnKRWq1mwYAFVVVWsWrXKbcctKSlh3bp1ZGdnM2zYMLcdt7+pItr+Dxy1dV6OxDdZ7E5e3VTImdlRTEv3/M8rwMSUMBaMjOH1zYWYbOKuQui+QZcoXK2tKIIC++VcI0aMYMqUKWzdupW1a9f2+Xl8aWkpb73zLi0uNX/aq2fGk2tZsafcTdH2L0VA2/+By2j0ciS+6btD1TSZ7PxyZmq/nveXM1NpsTj45kC/diQWBrhBt45CttmQNJp+O9/8+fOx2Wxs2LCB6upqzj///B/XQHSXy+Vi27ZtfLdqFQanmm+smVhRdXtxni86+X8g/6Sjn/A/647UEO6vYWpa/9xNnDQxOZSoQC1rjtSwZPzAnCgh9L/BlyisVhSa/luJqlAoWLRoEVFRUaxevZrnnnuOKVOmMHHixC6nttrtdg4dOsSmTZuora2lVhHGamsy1p/8t3RncZ4vUmhPJgqrlyPxTXtLm5iYEopS0fM6Xr0pCXOSQiExKTWMfaVNPT6vMHQNukShCAxEEeyeeejdJUkS06ZNIzMzkzVr1rBx40Y2btxIfHz8j61Q/f39USgUWCwWmpqaKC8vp7CwEJvNRkREBJdeeinnv1OIzOlvHF0tzvNJJ8q7y91YvT7UuFwypQ0mFozqebvY3paE+an0yAC+yq3E4XShUg66p8+CBwy6RDFstfsGlnsqPDycyy67jIaGBnJzczl69Cjbt29vd1bUyVaoI0eOJCUlBYVCQVxINeXtJIWuFuf5opOPnBT9+BhwoDDZnThcMq5ejGn1tiTMTwXpVMgyGG1OgvUiUQhdG3SJonXzZhR+fvjl5HgthrCwMM444wzOOOMMnE4nBoMBs9mMy+VCq9USFBSEVnv647H75medcrUI3Vuc54tc1rZHTlI7r3Ooszva7rKO1/R8oL+3JWFOOb+zLUGplQOnfL3gXYMuUVT/9XG0mZleTRQ/pVQqCQ0N7VY58pNXhL19/uxLZHPbG5ek1Xk5Et8T6q8hxE9NVFDPk2hciL7Pd51VzWYCtKrTKgYIQkcGXaJQhYfjqB+4c/cX58QPyMTwc44TtapU4WFejsQ3ZUUHcqC8ucf7ueOu80CFgczogAHVEEvwrkH3gFIZEYGzTqwG9jbniRXZyn5Y/DgQTUsPJ7e8mZoWS4/2601JmJ9qMNrYW9rE9HTPVy8QBo/Bd0cREYGxthZZlsUVkxc5autAklCFiTuK9pw/JpZ/rT7KxzvLuGNuz1bg9+Wu89NdZThdMueNie3V/sLQNOjuKDSJCbiMRpzd6BUheI6tpAR1XBySatBdi7jFsKhAZmVE8PqmQlos9n45p9Hq4JWNBUxNC2N4bP9OIRcGtsGXKE70oRC9ELzLVljo0Z4gg8F987NoMNl4+tu8fjnfM9/lU9Ni5b752f1yPmHwGHyJIrWtdo6ti/LfgufIsoytqOjH/wuhfWMSQrhheipv/1DMl/srPHqubw5U8frmQq6blsyE5K5n4AnCTw26RKGOj0fh74/l0GFvhzJk2YuLcRmNaLMyvR2Kz3vg3CwmJIfy26X7WJ9f65FzbD5Wx68/3MPYxBAeXjjcI+cQBrdBlygkpRLdqFGY9+/3dihD1snvvX7MWC9H4vu0KiWvXjeRYZEB/OrNHXywvcStXQE/3lnKDW/sICXcn9d/MRHdEF47YbXV0dCwmcrK5VRUfEx1zVcYDLm4XKJwZVcG5Uijfsxo6t98C5fVikKsDO535v25SH5+aIelezuUASHUX8MHN0/lzvd389CyXNbn1fKnRSOIDe596ZZqg4VHvzjEytxKZgwL58WrJhDs59nmSL7IZqujvPxDqmu+xGg82u42CoWOsLCZxMVeSkTEPCRp0F0/99ngTBTjx8Mrr2Leswf/qVO9Hc6QY9q2Df3YMR5vRzuYBOvVvHnDZF7eUMC/Vuez7ukarpycxHXTkkmLDOj2cYrqjLyztZj3thXjktsGzG+ZnTbkiv+5XDaKil+iuPi/uFwWQkKmMCz9AQKDRqPTxiBJGhwOAyZzIU1N26mt+Y79dasJCBhOdtajBAeP9/ZL8CmSO29z3WnixInyzp07e7Wvy2gkb+o0wn9xHVG/+52bIxM6Y6+u4dgZZxD1u3sJv/FGb4czIJU2mHh2zVFW7CnH4ZIZFR/E7IxIxiaGkBbhT0SAFq1agdXuot5opaDWyL6yJjYerWN/WTMqhcSisXH85qxMksL9vP1y+p3FWkXu/tswtOwnKmohaan34O+f1uk+LpeD6povOX7871it1aSn/47kpJsH3FosSZJ2ybI80e3HHYyJAqD4ul/gNBhIW7HcjVEJXWlatpzKhx8mdcVydNliGmZf1BgsrNhbzncHq9lT2oTT1fHvqlIhMTYhmHNGxnBRTjzRQUOzxpbFUsGu3VditzcyYvjTREXN79H+Dkcrh488TE3NSpKSbmRY+oMDKll4KlEMykdPAAFnzKbm6b9jKytDkyA6efWXlrVrUEVFoc0aeBVvfU1UkI6bZ6dz8+x0zDYnR6oMlDSYqG+1YXe6UCsVhAdoSAzzIzsmED/NoP117han08K+/TdhtzcxPuddgoLG9PgYKlUAo0Y+S746jJKSV9HrkkhIuNoD0Q4sg/YnK3D+Amqe/juGr78m4qabvB3OkOBsacG4fgOhV105oK7CBgK9RklOUig5SWINREcKCv9Ja+sRxo59rVdJ4iRJksjM/CNmSwlHj/2VsLAZ+PmluC/QAWjQjnBpEuLRjx2L4auvvR3KkNGyeg2y3U7QwoXeDkUYYszmMkpL3yQu9jIiwuf0+XiSpGB49pNIkprjBc/0PcABbtAmCoCg887DevgwlsNi8V1/aF62DHViIroxvb+aE4TeKCt/B4DU1LvddkytNoqEhGupqfkai8WzK+d93aBOFMEXLkLS6Wj84ENvhzLoWY8exbRjByGXXSoeOwn9SpZlaqq/Ijx8Djqde6vixsVeBsjU1Hzj1uMONIN2jAJAGRxM0HkLaf7yS6Lu+x3KwEBvhzRoNX74EZJaTcjFF3e6nSzLNDU1UVNTQ0NDA0ajEZvNhizLKJVK/Pz8CA4OJiIigqioKNTqobdITOgZq7USi7WCpOSbu9x2xZ7yHnWQ9PNLRq9Poql5J0n80p1hDyiDOlEAhF51Fc2fLqPpo4/EvH4PcTQ20rR8OYHnLmi3/4Tdbufo0aMcPnyYgoICjMb/9YpWKBRoNBokScLhcGC320/5WkJCAhkZGYwcOZIw0dtCaIfJVAhAgH9Gp9ut2FN+SnfA8iYzDy3LBeg0Wfj7Z/54jqFq0CcK/ciR+M+YQf0bbxJ69dUo9L0viyC0r+Htt5FNptNmlzU3N7N161b27NmDxWLBz8+P9PR0kpKSiImJISwsDD8/v1MeVdntdpqamqitraW8vJyCggLWrFnDmjVrSE1NZdq0aWRkZIjHW8KPWo35Jz7q/En609/mndJCFsBsd/L0t3mdJgq1KohWR2tfwxzQBn2iAIi47VaKr7mWpo8/Juy667wdzqDiNBhofOddAs85B21G2xWdxWJh/fr1bN++HZfLxciRI8nJySElJQVlF2U91Go1kZGRREZGMmLECACamprYv38/O3fu5P333yc+Pp758+eTlJTk8dcnDBwul7nTr1c0tf/1jj5/ktNlRVJoeh3XYDAkEoXfxIn4TZ5M3X9fIviii8RYhRvVv/wyLqORiNtuBSA/P5/PP/+c1tZWcnJyOOOMMwgJCenTOUJCQpg9ezYzZsxg3759rFu3jtdff51JkyZx9tlno9EM7V/ioS4sdDoAdntzp9vFhegpbycpxIV0/pTBbC5x+yD5QDOoZz39VNT99+NsaKD+pZe8HcqgYSsro+Gttwm+8EI0WVmsWrWK999/Hz8/P2666SYuvPDCPieJn1IqlYwfP54777yTqVOnsmPHDl555RXq6+vddg5h4PHzS0Op9KOpeUen2903Pwv9z8qs69VK7pvfcRUBh6OF1tYjBAUO7SnfQyZR6EeNJPjCC2l4621sJSXeDmdQqHn676BSEXrXnXz88cds3ryZiRMnctNNNxEf3/Ez377SarUsWLCAa6+9ltbWVl599VXKyso8dj7BtykUasLDzqCm5htcLmuH2y3OieeJJaOJD9EjAfEhep5YMrrT8Ymamm+QZTsREfM8EPnAMWiLArbHXl1NwcLz0I0eTdIbr4sB0T5oWbuWstvvIOzXd7MmOJj8/Hzmz5/PtGnT+jWOhoYG3nnnHYxGI7/4xS88mqAE39XQsJk9e68jK/MvJCRc45Zjulx2tm0/D0lSMmXyVwPi/cJTRQGHzB0FgDo6mqj7fodp61aaly3zdjgDlrOlhao//wVNZiZbY2PJz8/nvPPO6/ckARAWFsb111+Pn58f7733Ho2Njf0eg+B9oaHTCQmZwvGCf2CxVrnlmCUlr2IyHSc97bcDIkl40pBKFAAhl12G38SJVD/5f9jKyr0dzoBU/fjfcNTVUXvjr9izdy+zZs1i0qRJXosnODiYa665BpfLxdKlS3E4HF6LRfAOSZIYnv04smwnN/d2HA5j1zt1or5+PQWF/yQqaiEREWe5KcqBa8glCkmhIPaJv4EsU3Hvvcg/WeAldK35s89oXrECzU03siY3l7S0NObOnevtsIiIiGDx4sVUVlayYcMGb4cjeIGfXyojRzyDwZDL3n3XY7PV9eo4NTXfsD/3Vvz9sxie/cSQv5uAIZgoADSJicQ+9ijmffuo+de/vB3OgGEtKKTyL4+inziB7dHRSJLEhRdeiELhGz9G2dnZjB07lk2bNomZUENUZOQ5jB71HC0tB9m2/QKqa76mu+OwdnsjR448Qu6BOwgIGMn4nLdRqbrfhnYw843fcC8IOvdcQi6/nIbXXqf5y5XeDsfnOZubKbv9dhQaDfLvfkdefj6zZs0iODjY26Gd4uyzz0alUrF27VpvhyJ4SVTUAiZO+BSNJpwDB+5k+45FlJW9i9lcflrScLlsNDXtJC//z2zecgblFR+RlPgrJox/D7Va9P44aUgsuOtIzO8fxnr8GJUPP9zWv2LcOG+H5JNku52y3/wGW3k5yW+8zmeHD6PX65kyZYq3QztNQEAAkyZNYvPmzZx55pmiPtQQFRg4nEkTV1BVtYLS0tfJy/8T8CfU6jC02hgUCjUOh+FE8rAhSRqioxaSnHJrlzWjhqIhnSgkjYaE556j6LLLKb39DpLffRdtWqq3w/IpsstF5R//hOmHrcT+7W84MjPJX7mSGTNm+OyK6ClTprBlyxb27NnDmWee6e1wBC9RKFTExV1CbOzFGI1HaWz8gVZjHjZrLS7Zjk4XT2TE2QQFjSMsbDoqlajY0BGPJQpJkqKBT2RZniVJkgooOPEH4C5ZlnM9de6eUIWGkvjSSxRfcw0lv/wlye++iyZBzMWHtpLg1Y//jebly4m44w5CllzEtm3bkGWZcT589xUUFERaWhoHDx5k3rx5YjByiJMkiYCATAICMr0dyoDlkTEKSZJCgbcA/xOfGgN8IMvynBN/fCJJnKRNSyXp9ddwmUyU3HAD9nIxbVaWZWqe/juN771H2A03EHHnHQAcO3aMsLAwIiIivBxh5zIzM2loaBDrKgTBDTw1mO0ELgcMJ/49FThfkqTtkiS9duIOw6fosrNJeuVlnE1NFF19Ddbjx70dktfITidVf/wjDa+/TuhVVxF1/31IkoQsy5SVlZGcnOztELt0MkZR2kMQ+s4jiUKWZYMsyz8t5bgDOEuW5cmAGljY3n6SJN0sSdJOSZJ21tbWeiK0TunHjiX53XeQHQ6Kr74G8969/R6Dt7ksFsrv+S1NH39C+K23EP2HR358dGM0GjGbzURHR7vtfCv2lDPjybWkPriSGU+uZcUe99zNRUREoFAoqKmpccvxBGEo66/psftlWa488fFOoN1pBbIsvyzL8kRZlidGRkb2U2in0mVlkfL+eygCAym+9jqalq/wShzeYK+spPiqq2lZtYqoBx8g6je/OeX5vsHQdoPorimxJzuOlTeZkflfxzF3JAuVSkVgYOCPMQuC0Hv9lSjekSRprCRJSmAxsK+fztsrmqQkUpZ+hH7iBCofeoiqx/+Gy2bzdlgeZdy+ncJLL8NWXEzCCy8Qfv31p21jsVgA0Ol0bjlnZx3H3EGr1WK1dlxNVBCE7umvRPEo8A6wF/hBluXV/XTeXlOFhpL0yiuE/eI6Gt95h6LLr8BaUND1jgOMbLdT889/UfKL61H6+5Py0YcEzmu/JIfJZAJwWy2l3nYc666T4yqCIPSNRxOFLMtzTvx9QJblMbIsj5Zl+feePKc7SSoV0Q89RMKLL+KoqqJwycU0vPUW8iApOmc5coSiK6+i/qWXCL54CanLPkU7bFiH2wcFBbn1/B11Fuuq41h32Ww2n13rIQgDyZAt4dETgfPmkvrZCvymTKb6iScpuuxyzLk+NcO3R1wmE9VPPU3hxZdgr6gg/tlnifvrX1H4+3e638lE0dzcecvJ7upNx7HucrlctLS0ECja3gpCn4lE0U3qqCgS//tf4v/1Txy1tRRddjkVDzwwoEqVy3Y7jR9+xPH5C2h4/XVCliwh/auVBM0/p1v7BwUFoVKpcNeMtN50HOuupqYmHA6Hz6/3EISBwOfWM/gySZIIWrAA/xkzqH/pJRreeZfmr74m9LLLCLvhBp9d0S3bbBi+/pq6//wXW1ER+pwc4p99Fr/xOT06jkKhIC4ujtLSUrfFtjgn3i2J4edOxhgXF+f2YwvCUCPuKHpBGRhI1O9+R/q33xCy+EIaP/qI4+ecQ/lvf4t5716fGUB1NDZS/9prHDv7HCoeeLCtttWLL5L8/ns9ThInpaWlUVFRQWtrq5ujda+jR4/i5+fn1jUfgjBUiTuKPlDHxBD72GNE3HEHDe+8Q9NHSzF89TWaYemEXHQRQeefj7qf36hkm43WLVtoXraclnXrwG7Hb+pUYv/6GP4zZ/a57tHw4cP5/vvvOXDgAFOnTnVT1O5ltVrJy8tj9OjRPtMrQxAGMpEo3EAdE0P0ffcRcdttGL7+muZly6l5+u/UPP13dKNGETBvLgEzZ6IbPhxJrXb7+e1VVZi2b6dl3TqMGzbiMhpRhoYSdtWVBC9Zgi6r74PDJ0VHRxMXF8eOHTuYPHmyT74R79u3D7vdTk5O7+6aBEE4lUgUbqQMCCD00ksJvfRSrIWFtKxaTevatdQ99zx1/34OSadDP3Ys+tGj0KQPQztsGJrkJBSBgd260pdtNuxVVViPH8d2/DiWI3mYd+/GXlHRdv7ICIIWLmxLTDNmIHloauj06dP55JNPOHjwIKNHj/bIOXrLbrezadMmEhISSEhI8HY4gjAoiEThIdrUVLQ330TEzTfhqKvDtHMnpl27Me/aRcNbb5/Sq1vS6VBFRqIMDUWh0SDpdKBUIFttyFYrLqMRR20tzp9VQlVFR6MfN46w63+BfvwEdCOGI/XDFf6IESOIjo5m9erVZGVl+dRahS1btmAwGLjoootEeXFBcBORKPqBKiKCoAULCFqwAADZ4cBWUor12FHs5RU4amraEkFTE7LVitNgAIcDSadDodehDAtDPz4HVWQk6uhoNGnpaIelo3TzArjuUigULFy4kDfeeIPVq1ezcGG7NR77XXV1NRs2bGDEiBGkpooGVILgLiJReIGkUqFNSx3Q3fSSk5OZMmUK27ZtIzk5mZEjR3o1HovFwscff4xOp/OZxCUIg4XvjUQKA8bZZ59NQkICy5cvp6SkxGtx2O12li5dSkNDA5dccgkBAQFei0UQBiNxRyF0m8luorC5kApjBU3WJmxOG5oJGpRrlbz9ztvMXzKfCdkTUEj9d/1hs9lYunQpBQUFLF68WDxyEgQPEIlC6JBLdrGrehdrS9aytXIrx5uOI3P6YkJ9iJ5ZVbP4fOnnPBn9JFHDopgZP5Nzks8h0s9zfUWam5v56KOPqKio4IILLvDpPt6CMJBJvrKK+OcmTpwo79y509thDElmh5lP8j/hvcPvUd5ajlapZUL0BMZFjiMzNJP4wHhCtaHoVDqcspNWWytl9WVs/nozxhojdaF1bAncglPlZG7iXH456peMiRzjtvhkWebw4cN88cUXOJ1OlixZQnZ2ttuOLwgDlSRJu2RZnuj244pEIZwkyzLfFn/L0zuepsZUw/io8VyWdRlzE+fip/brcn+n08n69evZuHEjOr0OMuBLy5c02Zo4O/lsHpz8IFF+UX2KsaamhtWrV5Ofn09sbCwXX3yxKPwH4LRDcxk4LKAPhYBoENODhxyRKIQOWRwW8hrzKG0ppdXWilKhJEIXwbCQYSQEJnRrPYHZYebPW/7MV4VfMSJ8BPdNvI+JMb37eauoqGDlypWUl5cTEhqCM9HJxy0fo9Kq+NvMvzE7YXaPjifLMsXFxWzbto3Dhw+j1WqZPXs2U6dORalUdn2AwcrlhEMrYM+7ULQZnD/p5ucXDpnnwuQbIU6sUB8qRKIQTuGSXawvXc+nRz/lh4ofsLnab9Ua6x/LgpQFXJ59OfEB7VdpbbY2c9vq2zhQd4Dbx93OjaNvRKXo2/CVLMscOXKETZs2UV5ejkKpoDmgmeOq41wy+RKunHhlp2/yLS0tlJeXU1BQQF5eHs3Nzeh0OiZNmsTUqVPx76J3xqBXfRCW3wJVuRCSDNnnQ/QIUPuBqR7KdsCRr8DWAmOvgnP/D3TeWXcj9B+RKIQf7a/dz1+3/pXDDYeJ8ovinORzmBgzkdSgVIK0QThcDmpMNRyuP8ym8k1sKt8EwBXZV3BXzl2nPEayOq3c+O2NHKw/yNNnPM2ZSWe6Pd6Kigr279/P4cOHf2x6JCkkwkLDCAgIQKvVIkkSDocDk8lEU1MTZnNbO1SVSkVqaiqjRo1i+PDhPrUK3GuOrYGPrgFNACx4AkYugfZW5FsMsOmfsPlZiMiE61ZAYEy/hyv0H5EoBGRZ5p1D7/CPXf8gQh/Br8f/mnNTz+3y6r/KWMXL+1/mk/xPSA5K5t/z/k1qcNs00ie3P8l7h9/j72f8nfkp8z0ef3VdNY9+8ygtDS2cGXEmLqsLm82GLMuoVCr8/PwICgoiIiKCuLg44uLiUHugkOKAVbEHXj8XwofBNZ90742/4Hv44CoIS4NffQuaIX43NoiJRCHw/J7neWn/S5yVdBaPzniUQE3P2nzuqNrB79b/DoDX57+OzWnj8i8v54rsK3h4ysOeCLld1cZqLvrsIsZEjuG/Z/+338474Nkt8J/p4LDCzd9DQA+mHh9dDe9dApNvgoVPeyxEwbs8lSjEyuwB4rNjn/HS/pdYkrGEZ+Y80+MkATApZhJvn/s2CknBHWvu4MV9LxKgCeCunLs8EHHHov2juXnMzWyu2MzBuoP9eu4Bbedr0HAcLnyuZ0kCIOOstiSx41WoO+qZ+IRBSyy4GwCqjFU8vu1xJsVM4g9T/9Cnlc/JQck8O/dZrvvqOipaK7h6+NW9Sjor9pTz9Ld5VDSZiQvRc9/8rB61NL0482Ke2/McXxZ8ycgI79aJ8jW1tbUcOHCAoqIiamtrMZlMKCSJX/MaZn0GRbV+jIpu7Xmpktn3w663YNtLcN7fPRO8MCiJRDEAvLD3BZwuJ4/NeKzPs5EAxkSO4eyUs/mm6Jsfxyp6YsWech5alovZ7gSgvMnMQ8tyAbqdLAI1gYyPHs/OavF48aS6ujq+++478vPzkSSJ2NhYsrOz8ff3J9BwlOB9zWxTns2Wb75h1apVTJgwgblz56LX67t3goBIyDoXDn/e9vhJrLMQukkkCh/XaGnky4IvuTTz0g6nt/bGyIiRfFP0Ddsqt3FZ1mU92vfpb/N+TBInme1Onv42r0d3Fekh6aw4tqJH5x6sdu3axVdffYVKpWLu3LlMmDDh1DuGH/JgH5xzy+OMMyvYunUrO3bs4NChQ1x22WXsrld27w4v7Yy2tReNRRAm6mIJ3SMShY/7vvR7HC4HSzKWuPW4a4rXALClfAsOl6NHdyoVTeYefb4jOqUO608XiQ1R69evZ926daSnp3PRRRe1/0ipsQi0wRAYQ1QgLFq0iEmTJvHxxx/zxptv8b09g3J72yPETu/wIrL+dzyRKIRuEoPZPm5X9S7CdGFkhXbd93rFnnJmPLmW1AdXMuPJtazYU97hthdnXgxAq6OVY03HehRTXEj7jzo6+nxH6i31hGhDerTPYLNv3z7WrVvH2LFjufrqqzsed9j1BlibT/lUbGwsN954I82ynumKfEIk049fO3mHdxrtifEom9FdL0EYAkSi8HGlLaWkBad1WYbj5LhBeZMZmf9dVXaULIaHDf/x42JDcY9ium9+Fnr1qauq9Wol983vOpn91P7a/d1KgINVc3MzK1euJDk5mUWLFqHorI1tB2U4/Pz8+MacjgMFs9SFSD+p7tvuHd7JBKHW9SV0YYgRicLH7a7Zzf7a/V1u19m4QXsyQjOI1LdNsTTae3Z1uTgnnieWjCY+RI8ExIfoeWLJ6B6NTxxpOEJBcwGzEmb16NyDyfr163E6nVx44YVd16zKPr/tb1PDaV8KDQlmmz2ZcIWJNEX9j59v9w6v4Xjb3yHJvQ1bGILEGMUA4JJdXW7T03EDhaRgfsp83j38LlXGqh7HtDgnvkeJ4ef+s/c/+Kv9OT/t/F4fYyAzm83s27ePnJwcwsLCut4h4cQaqsINMHLxKV+6b34WDy2z0uCqYLiqmuO2iI7v8Ao3tlWXDUvv+4sQhgyRKHzc2clnc6ThSJfbxYXoKW8nKXQ2bpATlcO7h9/ls+Ofcf3I67tVStwdvin6hrWla7k7526CtcE/ft7caqPscCPVxQaaa8xYTXaQQaNXERiuIzIpkMThYQSGDfzHJkePHsXpdHbZbOnkepXqpla26UKwb3iDmJ8lipMJe+nKGjIdhaQHS9y1oJ07PEszHPkSRixuvzaUIHRA/LT4uFERoyhtKaXaWN3pdr0ZNzhQfwAFCipbK3lw44PYXXa3xNyZ/bX7+ePmPzImcgzXj7oeWZYpO9LAly/s4437N/Pdawc5sL4cQ50ZhVJCoVJgbLaSv62Kde8c4e2Ht7Ds6V0c3VmNy+Wb5We6o6ysDI1GQ1xcXIfb/HTcyYGSt+1nElP9PevXfHXatotz4vn7L88G4NkLktq/29vyHNha21ZoC0IPiDsKHzcncQ7/3PVPvir8ihtG3dDhdiffGLq7WtrhcvBN4TdMjZvKGQln8MT2J/j12l/z1OynCND0cMVvN20u38y96+8lQh/Bs3OfxVTvYP0HByk91IA+SEPO2UmkjYskMikAhfLUaxjZJdNQZaRofx2Ht1Ty3asHiUgsZs5V2USnDrzy2U1NTYSGhnY6gP3zcafXnOdypWotSZvug5mz/jeD6YTw8PAfj32asp2w6V8w+jKIG+eGVyAMJSJR+Li04DQmxUzi3UPvcnnW5Z0+HurJuMFXhV9RaazkgUkPcGbymagUKv627W9c9uVlPDbjMSZET3DXS8DsMPOfff/hzQNvkhGawYtnvojhiMzqN7eDBDMvzWDk7DhU6o4HdCWFRHhcAOFxAYw/J5mju6rZ8ulxPn16F9OXpDP2zMRuNWjyFfn5+V1u8/PxpVb8uMd+O++on4AProQr3j+lx8TJKrsOh+PUA1XlwvuXQ1BcW18KQegh8ehpALgr5y5qzDU8t+c5txyv3lzPMzufYVT4KOYmzQXgsqzLeG3+a7hkF9d/cz2//f635DW0P2Oqu6xOKx/nf8yiFYt448AbLMlYwjvnvkP9Hidfv5xLaKw/V/5xCmPPTOw0SfycpJDInBTDlX+aQurYCDZ/cowtnx7DVyshtychIaHL3hrtjS/94BrJ45q7oXgLvDKv7e8TLBYLAFqttu0TTgdsfwVeOwdUWrh2Ofh1Y+BcEH5G3FEMADlROVyZfSXvHn6XEeEjuCD9gl4fy+Kw8Nvvf0urrZU/n/3nUwoMToiewLJFy3jj4Bu8ffBtVhWvYmzkWM5NPZcZcTNIDkru8qrd7DCzp2YPa0vW8m3RtzRZmxgVPoonZj7BxJiJFOytZd17R0gaEca5t4xGpel9K1OtXsWCm0axcelR9q4uRR+oYfz8gTHtMzk5mYqKCux2e4f9NtpmM+We8vhJr1YyduHNEHIWLL8V3jgXEqfA8AswyGEkUk6ScR+s+gYOfArNpZA2Bxb/F4Ji++nVCYON6EcxQNiddm5dfSs7q3fy+ym/59LMS3v8qKXJ0sQ939/DrupdPDX7KRakLuhw22ZrMyuOreCz459xtLGtLHWoNpSM0AziA+IJ0YWgU+pwyS6MdiPVpmqKDEUUNBXglJ3olDrmJM7h0sxLmRQzCUmSaGmw8OGj2wiJ9uOie8f3KUn8lOySWfX6QY7uquGie8cTNyzELcf1pPz8fN5//32uvvpqMjIyOtyu0yq9NlPbiu3d70Dt4VN3VKggdTZMvhkyF/hEAUC7vRmzpRS7rRFZdqBQaNFoI9HrklAqtd4Ob1AQjYsETHYT966/l03lmzg7+Wzun3Q/Mf5ddziTZZk1JWt4YvsTNFoa+euMv7IwbWG3z1tiKGFb1TZya3M53nycytZKmqxNP86S8lf7E6mPJDEwkeywbHKicpgYMxG96tRHJ1/9Zz+lRxq54pHJBEf2rNxHV2wWBx/9dTtKlYIr/jD5tMFwX+NwOHjmmWdISkriyiuv7Pvxmir45KUniQwN5MzzLobILK93snO57NTXr6e29jsam7ZisbRfJUCSlAT4ZxMWNoOoqIUEBo4aUONNvsRTiUI8ehpA/NR+PD/ved44+Ab/2fsfvi/9ngUpC1iQuoCJ0RNPGeiWZZmyljI2lm9k2dFl5DXmkRGawXPznmNE+IgenTcpKImkoCQuzbz0lM+7ZBcSUrd+qauLDBTuq2PKojS3JwkAjU7FjEsy+Pq/ueRvryZ7mvceszicLiqbLRgsdvRqJfGherSqU++eVCoVU6ZM4fvvv6e4uJjk5L49Mtt6oIAj5jAmX3odxKf16Vh95XSaKSt7h5LSN7DZalCpggkNnUZC/NXo/VLQqMORFGqcThM2ay1G41GamndRUvo6xSUvExQ4hpSUO4iIOFMkDB8hEsUAo1QouXH0jZyXeh6vH3idLwq+4IuCL5CQiPaPJkgThMPloNZUS4u9BYCs0Cwenf4oF6Rf4JZ+Fif1pIHSge/L0OiUjJmX4Lbz/1zq2AjC4vzJ/b6s3xOFyyWz6nA1S3eUsuV4/SnjCiqFxPjkUJbkxHPR+Pgfk8a0adPYs2cPy5cv5+abb8bPr3cLHsvLy1m3bh3Z2dmkpXk3SdTXb+RI3u+xWMoJC51JYvbjhIXNQqHouu+53d5MVfXnlJa+zv7cWwgLnUF29t/Q6z33MyN0j3j0NMBZnVZ2Ve1iX90+ylrKaLG1oFKoCNeFMyxkGJNiJ5EW7O0rTBev37uRtPFRnHnd8K536IM9q0rY8ukxrv3rNIIi3H/n0p4jVQbu/2Q/+8uaiQ3WcfaIaEbGBRGs12CyOcirbmHt4RqO1rSSENpWF2tWRludrdLSUt58802io6O55pprepwsKisreeedd1Cr1dx88834+3vncZMsyxQU/ouioufx8xtGdtajhIZO6dWxXC4H5RUfcPz435EkBaNGPkt4+Gw3Rzw4iTEKYcCqLjLwyZM7mX/TKIZNiOrWPr1ttVpf0cqHj25n3nXDGT7d83cVqw5Vc9cHuwnQqnh44XAuHBePUnH64xJZltl4tI6/fHGQ47VGfr9wODfNbkvgeXl5LF26lMDAQBYvXkxKSkqX53W5XOzZs4dvvvkGvV7PddddR0REhLtfXrfIssyRvEeoqPiQ2NhLyMp81C2D02ZzCftzb8NoPMaokc8SFdXx5AuhjRijEAashopWACISu7fiuy+tVkNj/FGqFT+e05M2H6vjtnd3MTI+mFevm0hkYMdvjpIkMTszki/vmsXvPt7H418dRqdWcO20FLKysrjhhhv45JNPePPNN8nMzGTixImkpqaeNnXWZDKRl5fHtm3bqKqqIiUlhSVLlhAU5L3V6QWF/6Ki4kOSk28jPe1et40r6PVJTBj/IXv3/ZIDB+9hvCaCkBC3vwcK3SASheBxJoMNgICQ7l1l9qXVqkIh4R+s+fGcnlLbYuWuD/aQHhnAO7+aTJCu62fwAHqNkmevGIfV4eQvXxxiXGIooxOCSUhI4Pbbb2fr1q1s3bqV/Px8lEolERER+Pv743K5MBgMNDS0lRmPiIhgyZIljBo1qvM+Fh7W0LCZoqLniY29xK1J4iSVKpCxY15lx84lHDhwN1OmfI1aHdz1joJbiUQheNzRnTUAKNXde0Pra6tVhVKBy+nZR6pPf3uEFoudj26e2u0kcZJKqeCZS8dx9j/X88hnB1hx+3QkSUKj0TB79mymT59OQUEBRUVF1NXVYTKZkCSJ6Ohoxo4dS3p6OvHx8V6fEeRy2TiS9wh+fqlkZf7FY/Go1cGMGvUsO3cu4XjBP8jO+otHziN0TCQKweOyJsewpewYVpMDnX/Xb6q9KZn+UzazA7XOPYv52lPRZObT3eVcNy2ZjOjArndoR7Cfmt+encmDy3LZfKyemRn/G19QqVRkZmaSmZnprpA9orJqOWZzCWPHvIpS6dnS70GBo4iLu4KKig9JSb4Fna7jqruC+/n2qiRhUAiOanuDb6oxdbFlm760WrWZHZgMNoLCPTfj6bO9FThdMr+ckdqn41w0Pp5gvZpPd5e5KbL+VV7+PgEB2YSHz+mX8yUn3YIsO6mo+Lhfzif8j0gUgsdFJbdddVcea+7W9n1ptVpZ0HaOqJTeXel3x8ajtQyPDSIxrOOprCv2lDPjybWkPriSGU+ubbd3uValZF52FBuP1g6ogoYAFksFLS0HiIlZ3G+PwPT6eEJCJlFb912/nE/4n249epIkaQaQA+iBUuA7WZZPb94rCO0ICNURHu9PwZ5acs5O6tY+vW21WrC7BpVW6dF6T4crDSwY1XHplJ7M2hqbEMzyPeXUtliJCho4nfuam/cAEBo6rUf79Xba80mhodMpLPwXDkcrKpVn+qYIp+v0jkKSpFmSJH0MjAI2AZ8DTcDfJEl6VJIkzz0IFgaVrKmxVBU0U1Ns8Ng5LK128ndUM2xClNsKDv6c3emi0WQnLrjjR1udzdr6uZgTx6lpsbo3UA8zmYsA8Pcb1u19ftqxT+Z/CbS9u62O+Pu3nc9sLu1JuEIfdZgoJElKBS4ArpBl+SVZlvfKspwny/I3sizfCnwL3NVfgQoD24iZcWj9VPyw/LjHHrPs+KoQp91Fzlndu2vpjUZj27Tb5Z28ufVk1pZG1fbYxjHA2rpazG3jKj0ZxO5JAu2IWtW2XsThaOn2PkLfdZgoZFkulGX5flmWnQCSJE3+2dc3y7L8Lw/HJwwSWr2KyRekUXakkUObKtx+/IqjjeSuK2PkrHjC4jxXxiIyUItGpeCsEdEdbtPR7Kz2Pt9gbKvAG6zv2RRbbzMY9gH0KOn3ddpz2/nauvcp3FizTOhaTwazb5Mk6QdJkh6SJCnRYxEJg9boM+JJHB7Kxo+OUnG00W3Hba418c3LBwiK1DNtSbrbjtseSZJIDfcnr6rjK9qezNo6Wt2CWikR382pv74iLv4KAGy2mu7v04ME2hGLpe0iQ6PpXikYwT26nShkWb4BmA0cBtZKkrRGkqRzPBaZMOhIColzfjWKoAgdXzy/n5KD9X0+Zn15K8uf2YPsgvNuH4NG5/krzUmpoewoasDys8coJ/Vk1tb6/FpyEkPRqAbWBMSAgLbijgbD/m7v05dpzycZWnJRqQLFOop+1u3fqhOPnq4BpgEfn/jzMiDmqgndpgtQc+E9OXzx7318+fw+JpybwsRzU7q9avsk2SVzeEslG5fmo9GruPCeHEJj+qdy6sLRsby7tYQv9lVw6cT2b667M2trf1kTR6pa+PMFPesP4guCg8agVAZQW7eayMizu7XPye9Hb2c9uVwO6urWERo6HakHJe6FvuvJ5dftwNvAr+UTDyYlSXrAI1EJg5p/sJYl941nw4f57PyqiPwd1UxYkEzmpOguZyu5XDLFuXXs/KqImuIW4rNCOOv6kQSE9l8rzWlp4WTHBPL8umNcMDYOnbrnM6xkWeaZ7/IJ0qlYMmHg9VtQKLRERZ1LdfWXDEt/AI0mrFv79XbaM0Bd3WpsthpiYi7s1f5C73VYZlySpDggQZbl7R18XQfMkmV5lScCE2XGh4aSQ/VsWXac+rJWNDolSaPCiU0PJjTaH12AGkkBVqODphoT1YUGig/UYzLYCAzTMen8VLKnxiC1U9bb0zYereXa17Zzw4wU/nTByB7v/9GOEh74NJc/nD+CX83s2wpvb2ltzWfb9oUkJl5PZsYjHj2Xy2Vn+44LkGUHU6d8i5iZ3z5vlBmvBe6XJOlM4HlZln8cvZMkaTZwN/BndwckDC1JI8JJHB5GeX4T+duqKDlYz7Gd7Q+Qav1UJGSHkTEpipQxESi92Bd7VkYk109P4Y3NRcSH6LlxVvebQ605XM0jKw4wc1gE109P8VyQHhYQkEl83BWUlr5FVOQCj5YALyx6DqPxKGPHvCKShBd02bhIkqSpwHWADpBO/NkBvCbLssVTgYk7iqFJlmWMTTYMdSYsrQ5kZDR6FUHheoLCdV65e+iIw+ni7g/38FVuFVdOTuThhcMJ7KSSrN3p4uUNBTzzXR4j44J598YpA25a7M85HC1s37EIp9PMxAkfo9e7f0JkdfVKDhy8m9jYSxgx/P/cfvzBRHS4EwQf5HTJPP1tHi9tOE6Yn4ZrpiazYFQMmdGBKBUSsixTXG9izZEa3vmhiKJ6E+eNieX/Lh5DgHZwrAVoNR5l167LUSn9GTfujR9XT7tDVdXnHDp8H0FB48gZ95bHq9QOdF5LFJIkjZFlef+Jj6WfDGRfKstyh2UcJUmKBj6RZXnWiX+/BowAVsqy/NeuAhOJQhhI9pU28c/V+azPr0WWQa2UCNSpMVodWB0uAMYmhnD3vGGcObzjxXoDVUvLQfbsvQGXy0p21mNER1/Qp2KBTqeF48efprTsTUJCpjB2zEuoVJ4r9DhYeDNRrJVleV5nH7ezTyjwARAly/J4SZKWAItkWb5ekqTXgSdkWT7a2XlFohAGoqpmC1uO13G0ppUWix29WklKhD9T08JJjxzcRewslgoOHLibZsMewkJnkJZ+L8FBY3t0DJfLQU3NSo4X/AOLpYyEhGvJGPYQCkX/zWobyHylZ3Z3LxGcwOXAZyf+PQdYeuLj74CZQKeJQhAGophgHUvGD7zpru6g08UxYcJHlJW9Q0Hhc+zcuYSgoBxios8nLGw2fn6p7d5luFw2DC251NWtpapqBVZrFQEB2eTkvEtYD6vTCp7RnUQRI0nSVbQlieifftzRDrIsG4Cf/lD4AyerqDUA49vbT5Kkm4GbAZKSPFfYTRAEz5AkJYmJ1xMbezEVlZ9QUbGU/KOPAaBUBuDnl4xaHYZCUuN0WbBaazCbi5FlO5KkJCx0BlmZfyYi4kyxqM6HdCdRfARktPPx0vY3b1crbb0sAALooHSILMsv07bam4kTJ/rEKLssOzEaj2M0HsVsLsVmr8PhaEF2OZAUKlTKADSaCPT6RPz9M/DzSxcFy4QhT6UKJCnxBpISb8BkKqKxaRstLYewmEuw25twyQ6UCi3+fqlERpxJUNAYQkOnolaHeDt0oR1dvqPJsuyOTua7aHvctBUYC3S/rrAXWK011NR+S3399zQ17cDpNP74NaUyAJUqAElSI8t2HI5WnM7WH7+uUOgJCZlIRPgcIqMWoNN23OCmu5yyzMbGFlbXG9jVbKLYYsXgcKKUJCLUKjL8dEwN8ee8yBAy/MWsEMG3+Pml4OeX4u0whD7w6PRYSZK+l2V5jiRJQcBGYA1wLjBVluVO+2L292C2LMs0NG6mrPQt6uq/B1zo9SmEhc0gOGgcAQHZ6PXJqFSn1xNyOk2YzCW0th7B0LyXhsYtmEzHAYnwsFkkJFxHePicHs8Cccoy71fW82xxNWUWO3qFRE6QP8P8tISolDhkqLHZOdRq5pCxbUnLzJAAHkqLZUJw/9Q9EgTBdwz4dRQnZkKdDWyQZbmqq+37M1E0NGzh+PGnMbTsR6OJJC72EqJjLiTAP6PrnTtgNBZQXf05FZUfnxicG8Gw9PsID5/drf2LzVZuO1TMboOJCUF+3JoYxVnhQeg7WI1cbbWztKqBl8tqqbU5uCkhgj+kx6FRiOe8gjBUDPhE0VP9kShstjry8v9CTc1X6LRxpKbeRUzMhW6diudy2amu/oLCwucwW0qIjDyHrMy/oNV2XE//QIuJy/cVYJddPJmZyEVRId2+GzE6nPytoJLXyuuYFRrAm6NT8VeKkgeCMBSIROFmDQ2bOXDwHhyOFlJSbic56WaUSs/N1Xa5bJSUvE5h0b9RKv0YMfxpIiLmnrZdqcXGgp356BQSS8elk+7XuzGHjyobuOdICWeFB/Hm6FQUfVj8JAjCwOCpRDEkn0uUlr7Fnr3Xo9GEMXnSCtJS7/JokgBQKDSkpNzK5ElfoNXGsG//TRQXv3RKK0mXLHPnoWJsLhcf9SFJAFweG8ZjGfF8V2/g1bJad7wEQRCGqCE3j7Og4FkKi/5NRMRZjBzxj3YHpz3J3z+diRM+5tDh+zl2/Cls9kaGpT+AJEl8XtPEtmYj/8hOZFgfksRJv4yPYE29gacKq7g0JoxQ9ZD77xb6U3MZHPkKyne2fey0g184RA2H9HmQPAPEmNmANKT+14qLX6aw6N/Exl7CmNEv9nuSOEmp1DNq5LPEx19NSckrFBY9B8DLZbUM89NyRUz3msB0RZIkHkmPo9Xp4sPKBrccUxBOU3cMProG/jkKvr4PCjfQVvBKD03FsOXf8Nb58FwO7P2g7WvCgDJkLjFrar/l2PH/IzrqfIZnP+H1VZ+SpCAr88+4nBYKC5/Fpk5gtyGN36fFdjmesGJPebfbSY4I0DMmUM9Xtc3cliQa0gtuJMuw8zX45iFQamHWvTDuKghPP3U7ayvkfwNbnoMVt0LuUrj4NfBzzwWR4HlD4o7CbC7l0KH7CQoay/DhT3k9SZwkSQqysx8nJGQyZUf/QKxczoyQzgvHrdhTzkPLcilvMiMD5U1mHlqWy4o95R3uMyMkgH0tJpziSk5wp7WPwcp7IW0O3LULzvzD6UkCQBsAoy+Bm9bBec9A0WZ47Rxoqe73kIXe8Y13TA+SZZnDRx4GYNTI5zw+aN1TCoWaUSP/hUvScDMvkKDr/Cbv6W/zMNudp3zObHfy9LcdL3ZP1GmwyTINdodbYhYEdrwKG5+BCdfDlR9BYDdKpysUMOlGuG4FGCrg/cvA7rHeZ4IbDfpEUVP7NY2NWxg27AH0+t41dfc0rTaarzU3kUkeDTWfdbptRZO5R58HUJ54lOUUNxSCO9Tmtz1uyjgHzvtHzweok6fDxa9C5V5Y97hHQhTca1AnCll2UlDwDwL8s4iPu9zb4XRqbvoVHGcYlcXP4XLZO9wuLkTfo88DVFrtKCUIF7OeBHdY9UdQ6eHCF0DRy8Wc2Qsh5xrY+iI0FLo3PsHtBnWiqK1bjclUSErqnT7fkH1MkD/LuAyXrYKamq863O6++Vno1ae+Fr1ayX3zszrcZ1tzK8P99ah9qN+0MEDV5kP+1zD9Tgjo4+SIuY+ApICt/3FPbILHDOpEUVH+IVptDJER53g7lC4l6zSY/KbRKMVSXvFhh9stzonniSWjiQ/RIwHxIXqeWDK6w1lPJWYrW5uMnBMR5KHIhSEl9+O2N/cJ1/f9WEGxkH0eHPgUXM6utxe8ZtA+i7Dbm2ho3EJS4q8GRH8ISZL4ZWIUa/Nms6RpKVZrTYf1oBbnxHeYGH7u70VVqCWJa2LD3RmuMFQVfA/xEzu9m+jJ9G2yFsLB5VB9AGJ71jZV6D+D9o6iofEHZNlBROSZ3g6l266KDaNWNwMJmeLaDX0+3je1zSytauSWxEjidBo3RCgMabIM1Qchvt0GlUAvpm/HnThW9UH3xyu4zaBNFM3Nu1EodAQFjvF2KN2mUSj4w6gzMOHH54XfU2/r/XTWzY0t3HaomLGBeu5N6XvzJEHAagC7EYI6vpvt8fTt4BPHaumy84DgRYM2URhb8/D3z0ChULv92Cv2lDPjybWkPriSGU+u7XSxW0+NCQrAzz+LAHshC3fls6vZ2PVOP+GSZV4tq+XKfQUk6NS8OyYNXQc9LAShR06+mR/9rsNNejx9W3niTtfZ8Uw/wfsG7TuIxVqBXp/o9uP2ZmV0T8UEppCpbsIhy5y/+yh3HirmYGvH6yQA7C6ZlbVNLNiVzyNHy5kdFsjn4zOI1Lg/UQpDVGhK299JUzvcpMfTty0nGl1qA/sQmOBpvj/K20t2ezNqdajbj9vZrXV3B5i7olGHITmb+X5aNv8squaN8lo+qW4kXa9lWkgA6X5aQtRKnCdaoR5sNbO5sZUmh5MknYYXhiexJDq0x61XBaFTKi0EJUBjUYeb3Dc/i4eW5Z7yO9Lp9O26/La/w1LdGKjgboM2UTidJo8ctzcro3tKUqiRZTuBKiV/HBbHXclRLK9uZFW9gZW1TTQ6Tk1USToN50QEcX5kCPPCglANkfUSLpMdy9FGbGWtOA02kGWUARrUsf7oMkNRBvtWuZZBIWFCW60mWYZ2LkROXix1e9ZT0ca2v+MneCpiwQ0GbaJwuSzY7Y1uP25ciJ7ydpJCZyuje8rlsp7SjjVUreKXCZH8MiESWZYxOJw0O5woJYlwtWrIjUE4mqwYVhVj2lvTVpdEpUAVrAFJwtLSiGx1ggS67DCCzk5GE9d5oUWhB7IWwqHPoHgzpMxsd5NuT9+WZcj9FBIm9X3xXj+RHS4cdWYc9RacrTZkqxPZJSMpJRQ6FYpADapwHapwPdIgumAbtInC3z8Dl8vm9uP2+Na6F2y2OjTq9tc9SJJEsFpF8BAtx2HcU0PT8qPILvCfFIPf+Cg08YFIyrZfSlmWcdSYMO2rxbi1kprn9hB0VjKBcxMH1S+u1wxfBN88CJv+1WGi6Laj30Ht4bZSID5KdsnYSgyYDzdgPd6EvcIIrq6LpklqBeqEQHQZIehHhqOO9k7vG3cZtO82el0iZnOx24/b41vrXjCZitDpE9x2vMGiZUMZzV8VokkNIuzSLFRhp3cBlCQJdbQ/wef4EzgznqbPj2NYVYyjwULoxRkiWfSVxg9m3tNW7+nIyraV1b1hM7YlnLB0GH2Ze2N0A5fJTuu2SozbqnA2WUEhoUkKJHB2POoYf1ThehRBGhQ6JZJCQnbKuMwOnAYbjloz9opWrEXNGL4rxvBdMepYf/ynxuI/PhpJPfCeAAzaRBEQkE19w3ocDqPbO9n1ZGV0T7lcVozGfBLir/HI8QcqU24dzV8Voh8TQdjlWUjdeNym8FMTenkWyjAdLWtLUYXpCDozqR+iHeSm3Aa5n8CK2+HG1RCR0bP9XS747M62YoC/+AJUvrMYVLY7aVlfRsuGcmSbE+2wEILnp6AbHoaikxYAkhoUOhWqUB3a5CCgrey6s8WGObcO444qmpYfw7CmhOCzkvGbGD2gLloGXmrrppCQSciyk6bmHd4OpUeam/ficlkJDhGDeyc5jXaalh9FnRhI2GXdSxInSZJE0NnJ+OVEYVhTjK2i1YORDhEqDVz2NijV8Ob5ULmv+/vaLbD8Fji4DM7+C6TO8lycPWQtaqbqn7sxrC5BlxFC9G/GE3njaPxyojpNEp1RBmoImB5H1N05RNw4GlWIlsZlR6n97z7sde6bAONpgzhRTEap9KO2tuPFQb6opvZbFAoNYaEzvB2Kz2jdUIbL7CDs4gwkVc9/ZCVJIuSCNBQ6FYbv3P84ckgKS227G1Ao4dWzYcPfwdbFTMOizfDK3LZWqPP+ANPv7p9YuyDLMi0byqh9aT8AETeNJvzaEahj3PckQpIkdMNCiLxtLKGXZWKvMVPz3B7Mh+rddg5PGrSPnpRKHZGR51BdvZLMjEdQKv28HVKXXC4r1dWfEx4+D5VKzNQBkJ0ujDuq2gYE+/CLq/BT4z8tjpY1JTiaLKhCTh/fEHooajjcvB5W3tPWFvWHF2DkRW2D3KEpbesuWmugYg8c+RLKd7WV/7jqY8j0jYrOsizTvLKQ1k3l6EeFE3ppJgqt594WJUnCf3w02rRg6t89TP07hwi9OAP/ib5dZmfQJgqA+PirqKpaQUXFUhITr/d2OF2qrFyO3d5IQvxV3g7FZ9hKW3CZHPiN63u1Ur9xkbSsKcGS10jAlFhPht3vZFmmocJIdaGB5lozNrMDSSnhF6ghPN6f2GEh6Pw9sEo/IBIufxeKf4DtL8G+D2Dna6dvFzMazn0Kcq5tGxD3EYZVxbRuKidgehzB56f127iBKkRH5M1jqH/3MI2fHkXSKPEbE9kv5+6NQZ0oQoInEBIymaLil4iNvdTtg9ru5HRaKSp+kaDAMYSGTvd2OD7DVt42pqBJbr+fxsmSKienK58sqQKclixUEXoUfirs5YNnnMJitJP7fRmHN1fS0tDWf1qhlNDoVbicMjZzW2FJSSGRmB3K6LkJJI8Kd/+q/eRpbX8cNqg5BIZycNrALxwih7clFB9jPlBHy9pS/CZGE3xBWr9XMlBolERcO5zaVw/Q+HE+6ig/tz7ucqdBnSgAhqXfz85dl1BU9DzDhj3g7XA6VFz8XyyWcoYP/z9ReuMnnM1WABQB7V8N96SkiiRJKIO1OFvcv76mv8kumQMbytn6WQE2s4PE4aFMPC+F+MwQgn6y2MtudVJb2kLJgXrytlWx8oX9xA4LZs7V2YTFeuBNSaWBuHFtf3yYs9VG47KjqBMCCF08zGu/c5JaSfg1w6l+djcNS/OIuiPnxzVBvmTQJ4rg4BziYi+juORVIiLPIiTY92YTGQy5FBW/SHT0IsJCp3k7HJ/SurHzYos9LqmikKDr9VI+zWZxsOq1gxTl1pOQHcrMSzMIj29/TEutVRI3LIS4YSFMuiCVI1sq+WHFcT7+2w7mXptN5mTffjbuKYY1JbgsDiIvzezVBAl3UgZqCFmUTsP7RzDtqsbfB/9PBu2sp5/KyHgYvS6BAwfuxmqt9XY4p7DZGsg9cBcaTSRZmX/0djg+J3hhW7E4l6n93hw9rVbqbLZ2eHcyEFjNDj775x6KDzYw87IMFv16XIdJ4ueUSgUjZ8Vz5R+nEJUSxKrXD5H7fZmHI/Y9zlYbxh1V+E+I8ZkV0/rREagTAjCsL0Xuxsrv/jYkEoVKFcjo0S/gcBjYt/9GHI4Wb4cEgMNhZP/+m7HZahg96nmPVLsd6E4+s7WVtv9/dt/8LPRq5Smf66ikiqPJgqvV7rPPgbvicrr45qVc6kpbOfeWUYydl9irRyb+wVoW3T2OlDERbPgwn+N7anoVjyzLOOrNmA/V07q1AsP3pRi+L6X1hwrMh+qx15mRZd970zPtqwWHTMCMOG+H8iNJkgicEY+z3oKtyODtcE4z6B89nRQYOIJRI//N/txb2bP3BsaNfRW1OsRr8TgcLezbfzPNhn2MHvU8wcHjvBaLL9OmBCFplJj316LPDjvt6z0pqWLeXweALmtgJuSdXxdTdqSReddlkzq2b4PDSrWC+TeNZPkze1j71mGikoMIbKckys/JLhnrsSZM+2qx5DXgau284ZDCX4UuMwz92Eh0maE+sRrZkteIKkrvcxcMuhHhoJQw5zWgTQv2djinGDKJAiAiYi6jRv2bAwd+w85dlzN2zH/x8+v/Ovhmcyn7c2/FaDzGyBHPEBU1v99jGCgktRK/8VEYd1QRdHYyqtDT38y6U1LFZXPSurkcTUoQ6kjfmZ7ZXU3VJnZ9VUTGpGiGT3fPlbBKrWT+jSP54LHtbPwon4W3ddw2WHbJmPbU0LKuFEedGUmnRJcVhjYtuK32UagW6cTqZdnixNFkwV5lxFZowJLXgGlPDcowHUFzEvGbEO3VAVt7WQv6URF9Pk53p2V3l0KrRB0XgL2Du2dvGlKJAiAqcj45494g98CdbN+xmOysx4iJWdRv56+p+YbDRx4GZMaOeZXwcN8pYeCrAuckYtpVTePyY0RcP7JXV6WG74pxNtsIuzzbAxF63vYvC1GqFcy8tId1lboQFKFnwoJktn1WQFVhMzGpp1/J2quNNHxyFHtpC+o4f8KuyEI/KqLjQWCNEmWQBm1SEEyORXa6MB9qoGVDGY3LjtL6QwWhl2Z6pfy7y+LAZXKgCu9bW4CeTMvuCXW4DmuJ7yWKITFG8XOhoVOZNHEFAf4ZHDx0D/v23YTJ5NnSDmZzGftzbyf3wB3o9UlMmrhcJIluUoVoCT4vFWt+I80rC3r83Lt1WyWtm8rxnxbrc7f03WFssnJsVw0jZ8XhF+T+Anpj5iag0avYv/b0gW3T3hpqnt+Ls8FM6OVZRN2Zg9+4qB7NFJKUCvxGRxB1+1jCrsrG2Wqn5oW9GHdUufNldMvJSRHyz6ZU91Rn07L7QtIq2/qp+Jghd0dxkl6fwPjxH1Ja9iaFhc+ydds5xMZeQnLSjW59HGU2l1JS8hrlFR8hSQrS035HUtKNKBQDd+aNN/hPicVRa6Z1cwXOVjuhi4eh0Hf+4ys7XBhWF9PyfRm6rFBCzkvrp2jd69iuGmSXzIiZnhl81ehUZEyKJu+HSuw2J2pN2+SAlk3lNH9ZgCY1iPCrhqMM7FuSkiQJvzGRaNNDaPjwCI2fHsXZaiNobj9W9D1xM+oytz+Lrrs81elStrt8sgz5kE0UAAqFiuSkG4mJXkRh0QtUVCylouIjwkJnEBOziIiIM3s14G23G6ivX0dV9WfU129AklTExlxEaurd6HSDq3REf5EkieDz01AEaDCsKsJ6rImAWfH450Sd1vLUZbJjOlBH64ZyHHVm/CfFEHJhutfny/dWyaEGQqL9CPXg4Gva2AgObiin8lgTSSPCMe6spvnLAvQjwwm7Mtut3zulv5qI60fR+Ek+hm+LUehUBEzrnxlIyhAtkkbZ57U0nup06Wiw+GQL3yGdKE7SaqPIzvoLqSl3Ul7xAZWVn3Do8P2AgsDAkQQHjSMgIBu9PgmtNhqVKhCFQo3LZcfhaMVmq8FkLqa19QgGwz4MhlzAhVYbQ0rK7cTHX4VO63uLaAYaSZIImpuILjOU5q8LMXxThOGbIpRhOlQhWlBIJxrHmEAGdaw/4TeMRJ91+mypgaSm2EDqmO4PvvZmkDXmxCO5mqIWYoI0NC4/inZYiNuTxEmSUiL0kkxcFgdNXxxHHeuPNsXzjwUlSUId54+1pG9TUD3R6VJ2uLBXtPpkgUCRKH5Cq40kLfVuUlPuwtCyn/r69TQ2bqWyahlOp7HL/RUKPUGBo0hJuZ2I8DMIChqHJA3Mq1hfpokPIPLG0dhrTFgON2Ara8FpsIHDhSpCj350BPrsMNQJAQO+HIrN7MDSaickpnsztXo7yKrRq/AP1tBca6JhaT1KfzXhV3kmSZwkKSXCLs+i+t97aFiaT8w945F+tibGE3SZoRi+K8bRaGl3Fl13eKLTpeVoI7LNhTbT96Zvi0TRDkmSCA4aS3DQWEi9G1l2YrFUYjYXY7PV43C2IrvsSAo1KqU/Gk0Een0SOl0ckuT5H/ShpKW+jpID+6guOEZzTRXmFgMupwu1Tot/SBjhCYnEjRlOfPZIVOrBN+5jMbatU9AHdG98oCe1r35OF6AhoMqIo9FC+HUjUPh5/vup0KkIvTiDuldyadlUQdDcRI+f029cFIZVxRi3VhJ8bu/HI93d6bL1h0oUAWp0GSFuO6a7iETRDZKkRK9PQC/6WPcLWZY5un0Lu7/6nPIjBwFQa3WExMahDwxCqVRis1ioOpZH3g8bQZbR+vmTPXMOky5YQnBUtJdfgfu0NrZVhK0pMjB8etfjW30ZZFVIMpEGG5rUYHTD++9xnS49BF12GK0bywicGefxuwpVmA79mEhat1QQMD3OJ8YErAVNWPMbCVqQ0qMOjv1FJArBpzRVVfL1i/+kIu8QIdGxzLziOtLGTyIiMRlJcfovkNVkouzwAfJ+2MiBtd9yYO23TF1yBZMvuhSFYuDf3Z0cwO7uo6e+DLIGmhxoXDKBs+L7/ZFdwKx46l5pwHywvtPeI+4SfE4y5oP1NC4/RvgvRnj1EaXL6qRx2TGUIVoC3LSY0t1EohB8RnXBMT55/A/Isotzbr2bkWec2eWbvdbPj/QJk0mfMJlZV/2C9W+/xual71JbXMh5v74fhXJgJwtdgBqVVomhtnvTLns7yOp0uAi2OHD6ta247m/a1GCUQRrMB+r6JVGowvUEn5tC8xcFtKwrJWheP07R/QnZJdP4aT6OejMRN45GofHNn1eRKASfYG4xsPypR9Ho9Vz6yOOExPR8GnFgWATn/+YBYtIzWP/u6wR/+Dazr77BA9H2H0mSiIgPoLqbheJ6O8haW9pChEqCGH+vlNeQFBLaYSFY8hqQZblfrvADpsdhL2vF8F0xCn3/TdE9SXbJNH12DPP+OoIXpqJLD+nX8/eESBSCT9j+2SeYmpq45sl/dZgkujvtc+IFS6gvL2Pnl8sZd855BEV6/grVkxKyQ9n1dRHmFhv6bix6680ga8neWuIUEn7DQnoVozvqHqnjAzDtrsHVau/z4r7ukCSJ0Isz2qbofnYcZ7ONoHOS+6VwocvqpPHjPMwH6gmck0DALPcNinuC742aCENS/tZNpE2YRFRK+6unT077LG8yI/O/aZ8r9rTf2GjqksuRXS6Obv/Bg1H3j/TxUcgy5G3zTMkLl0umZGc1ALpe9Gfo6f9NR1QhbYPKJ7sa9gdJpSD8muH4T46h5ftSal/JxVHft9XVXbEWG6h5bg/mg/UEn5dG0PwUn5/GLRKF4HV2qwVDbQ11JUUdbtPT2jrBUdGodXoMdb3rteBLIhICiB0WzL41pTjtLrcf/9iuauxNbe1hZUfPj++uukfSiefzvYmhLySlgtAlGYRemom9opWqf+6i+ZsiXKbOS6j3lKPBQsPSPGr/sw/Z7iLixtFemTjQG+LRk+B1SlXbfH2tX8fVRHs67dPc2oLdYkYfGNT3AH3ApPNT+fxfe9mzqoSJC1Pcdly71cnW5QVEROrA4sDVi4J07qp7JJ9Igt6aHuo/IRrdsBCavi6k5ftSWreU4zc+Gv8J0b1evCk7ZawFTRh3VGE+UAcKiYDZCQSdmYhCO3DefgdOpMKgpVAqGTF7Hke3/4CxqRH/kNNXpvZ02uf+VV8DkDrO93qk90Zidhjp46PYsbKQxOFhRKf2PQHKssyGD/NoabRw5p1jkd85BL1ow+muukdOQ9tdjcIDFXK7SxmsJfyKbOxzEmnZUIZxZxXGrZUoQ7ToMkPRJAWijg1AFa5DoTv17VOWZVwmB446M/aKVqxFBqxHG3GZHEh6FQEz4gmYGY/KB9Zt9JRIFIJPmHLR5eT/sIkv/vkkFz/0F9S6U0sr9GTaZ3HuXn745H0yJk8nOm2Yx2PvL3OuyqKmyMBX/9nPRb8bT0hU3xow7f62mCM/VDFxYQrxI8OpCFRjr+y6VM3Puavukb2yFUnb1svC29Qx/oRdlkXI+WmYD9ZjPtyAaX8txu3/GyeS1AoknQpJKSE7XLgsDnD8L9EqAjXossLQjwxHlxXmk1Vhu0skCsEnhMXFM//23/DVv//OB3/4HQvu+O0pA9vdmfbpcjrZ880XbHjvTcLi4jn7lrv6/XV4ki5AzXl3jmHFM3tY/vfdLLx9DNEpPb+zcLlktn12nN3flpAxKZrJ57eVsdCmBGM93tTj6anuqHskyzLW481trW996Jm9wk+N/6QY/CfFILvktruFaiPOBgvOFjuy1YnsdCGp2pKGMkiDKkyHOta/rVKtD72WvpB8sfk5wMSJE+WdO3d6OwyhnxXu2cnXL/4TS0sLmVNnMHrefBJGjEKp6viaxtjUSP7WTez55gsaKytImzCZc2//LbqA/u+g1h/qK1r58vl9mJptTDovlXFnJ6LqZtmLhgoj379/hMpjzYycFcfsKzJRnBgTMO6pofGjPCJvHtPvDZ5s5a3UPLeHkMXDCJgqSvH3liRJu2RZnuj244pEIfgac2sLOz77hP2rv8FqMqLW6YlJG0ZobDz6oGAUSgV2q5WW+jrqSoqoLysBIDptGFOXXEH6xCmD5kquI5ZWO9+/n8fx3TX4B2sYNSeBjInRBEeePi7gdLioyG/i0JYKju+qQaNXMfPSDLKmxpzyfXLZnFT+bTu69GDCrx3Rny+Hho/yMB+oI/bhKV02pBI6JhKFMOTYbVaK9u6iOHcfNYXHaKquwtLSgiy7UGm0+IeGEh6fSFzmcNImTCYyKcXbIfe7siMN7Py6mPK8RgD8Q7SExvih9VPhcsoYm6w0VBhx2F1o/VSMmBFHzvykDqvRNq8qpmVNCZG3j23red0PbJVGav69m4AZ8YScPzC7EPoKkSgE4YT+KvEwkBjqzBTl1lNd2ExzrRmb2YFCKeEXpCEsLoD4rFASs0NRdVFLyGV1UPXMLhR6FdF35nh8AFZ2uqj5zz6cjRaifzsRpf/gKxXfnzyVKMQ9njDgiCRxuqAIPWPmJsDcvpXCV2jb+kPUv3GQxhXHCL0kw2Pfb1mWafqyAHtZK2FXZ4sk4cNEohCGNIfNdqIZkhONXo8uIHDIJyJ9VhiBZybRsqYEhZ+K4IWpbv+eyLKMYXUJxh8qCZgVj9/oSLceX3CvfksUkiSpgIITfwDukmU5t7/OLwgAFmMrR7dtoXDvTiqP5dNaX3fK11UaLRFJySSOGE3GlOnEpGcOycQRdFYSLpOd1o3luIx2Qi8a5raGQrLDRdMXxzFuq8JvQnSfuswJ/aM/7yjGAB/IsvxAP55TEAAwGZrZuuxDctd8h8NmJTA8koTskYTHJ+IXEoJCocRmNtFcW0N1wTF2rfyMHZ9/SlRqOjMuu4a08ZO8/RL6lSRJhCxKR+mvxrC6BFt5K6EXZ/R5gNtW3krjp/nYK4wEnpHQVhCvH6q1Cn3Tb4PZkiTdDtwBGIFc4BZZlh0dbS8GswV3Obp9C9+99BxWk5ERs+cx7pzziE4b1umdgsXYSt6Wjez8chlNVZVkTpnB2bfchc5/cK7N6Iwlv5HGT/JxGmzoR0cQeEYCmoTAHh3DVmmkdX0ppn21KPzVhC7JQD8i3EMRD10DftaTJEmTgDJZlislSXob+ESW5c9/ts3NwM0ASUlJE4qLi/slNmHw2vnFMta/+zox6RksuP0ewhN61snM6XCw88vlbFn6LqGx8VzyyF8JCO3/DnDe5rI6aPm+jNYtFchWJ+pYf3QjwtGmBaOO8T9tINplsmOvMmItNGA+VI+9vBVJrcB/WhxBcxPFWgkPGQyJQivLsvXEx3cDalmWn+loe3FHIfTVwfVr+ObFf5I5bRbn3vFbVOrez6opPbif5f/3KCExsVz52NOotbqudxqEXGYHpj01mPbVYisxwIm3D0lzou4R4LI4kW0n6j5JoE4IxG9sJP7jo1D4iZlNnjQYpse+I0nS48ABYDHwt348tzDEGGprWP3aiySOHMPCO+/ttARIdySOHMMF9zzIsif/zOaP3mHOdTe5KdKBRaFXETA9joDpcbhMdmylLdhrTDibbW1F8QCFVokyWIsqyg9NYqCY9joI9GeieBR4H5CAz2VZXt2P5xaGmK3LPgSXzILbftPnJHFSas5ERp85nz3ffMmE8y4iMDzCLccdqBR+anRZYeiyht6juKGm3+reyrJ8QJblMbIsj5Zl+ff9dV5h6LHbrBzZvIHhs+Z02S97xZ5yZjy5ltQHVzLjybVdtu+csvhSXE4nhzauc2fIguDTBm6BdEHoQE3BcexWC2kTpnS6XW96PQdHxRCRlELZ4QNujloQfJdIFMKgU1PctqYzKKLz1b697fUcGhNHS11t34IUhAFEJAph0KktaksUVlPn3dp62+vZ5XIhKcSvjjB0iJ92YdAZc+YCAEzNzZ1u11FP5656PTeUlxAcFd274ARhABKJQhh0IlPS0PkHkL9tc6fb3Tc/C/3P6hd11eu5tqSIxsoKkkaNc0eogjAgiEQhDDpKlYpR887h6NbN1JYUdbjd4px4nlgymvgQPRIQH6LniSWjO+31/MPH76PW6hg+8wz3By4IPkqsoxcGpckXXsLB71fz1XN/58pHn0Kj92t3u8U58Z0mhp/KXfcdR7dvYeYV16EP7J/ub4LgC8QdhTAo6QODWHjX76gvK2HZk3/G3NrSp+Md2rCWVS8/T9LocUxadLGbohSEgUEkCmHQShk7nvPuvp+qY/m888DdFO3d1eNjmFtb+O7l5/j6hX+QMHwUF977MAqle/oyCMJAIR49Cd3WXFPN8Z1bKTt8kPqyEoxNjTjtdlQ6HYHhEUQlp5I4cgzpE6f4TDnurGkzCYqM5Ovn/8GnT/yJhBGjGHfOeaTlTEKta7+wnyzL1JeVcHD9GnLXfIvNbGbSoouZcfm1bisHIggDSb9Vj+0pUT3Wd1TkH2Hrpx9QeOKKPDg6hqjkNALCwlGq1dgtFpprq6k+fhRziwGVWsPw2XOZuuRygiI6L6HRXxx2O/tXfcXOL1fQUl+LUqUiKjWdsPhE/IJDUCqVWM0mDLU1VBcep7W+DkmhIGPSNKZeciWRSSnefgmC0KUBX2a8p0Si8D6b2cS6t17hwLpV6IOCGXfOeYyYNZeQmNh2t5ddLqqOH+XAulUcXL8aSVIw4/JrmHDeYp9ZoOZyOSk7dJDCvTupOpZPY1UFZoMBl8uJRqcjMDySiKQUEkeMZtikqfiHhHo7ZEHoNpEohH7VXFPNsif/TGNFORMXLWHqksvR6DpfiPZThroa1r7xMsd3bmXYpKksvPs+1BqtByMWBMFTicI3LvMEn2Koq+WjPz+IqamRSx75K7Ovur5HSQIgKCKKC3/3e+b+4iaO7dzGiqcew2G3eyhiQRA8SSQK4RQOm40VTz+G1WTk0j/+jaRRY3p9LEmSGL/wQhbc9htKcvey7s2X3BipIAj9RSQK4RRbPn6P2qICzrv7PqJS0txyzJFnnMmkCy9h/+pvKNi9wy3HFASh/4hE4SGyLGMy2KgtaaG60EBzrQmX0+XtsDrVXFPFrpWfMWru2aSNn+TWY8+47GrC4hP5/u1XcDmdXe8gCILPEJPC3aym2MDBTRUU76/D2Gw75WsqjYL4rFCypsSQnhOJQulbeXrXV58BMP2yq91+bKVKzczLr+Xzf/yNYzu3kjllhtvPIQiCZ4hE4SbGZisbP8zn+J5aVFolKaPCiUkLJjBch0IhYWqxUVfSQlFuPd+9epDQGD/mXJ1NXEaIt0MHwOV0cnjTejImTyMwzDO9oNMnTSEgLJxDG9aJRCEIA4hIFG5QW9LCl8/vw2p2MPmCVMbMS0Srb+dbOwNmXS5TsLeWLcuOseIfu5l5WQZj5ib2f9A/U1N4HEuLgWGTpnZr+xV7ynn62zwqmszEhei5b35Wl8X1FAol6ROmcGjDWlxOpyiFIQgDhG89+xiAGquMfPavPShUEpc+OJFJ56W2nyROkBQS6eOjuPyRyaSMiWDjR0c5uLHjHs39peZEV7jYjOwut+1Nr+mTYjOysFstNFVX9jVkQRD6iUgUfeB0uPj2lQMolBIX/XY84fHdr2+k0alYcPMokkaGs+HDfOrK+lbdtK9aGuqQJAWBEV0/duptr2mA4Mi2znAtdXW9C1QQhH4nEkUfHNxYTn25kbnXZBMU0bMFaQAKpYKzbxiB1k/Fpo+PeSDC7ju6bQuy7EKh6PpxUG97TQOoNBoAHHZbF1sKguArRKLoJVmW2bemlNhhwaSOjez1cXQBasbPT6Y8r5G6slY3RtgzCcNHAeB0OLrctre9pgFsFgsAam37lVsFQfA9IlH0UkOFEUOdhexp7RfIg7Zn+TOeXEvqgyuZ8eTaDp/hZ02NAQkK99V6KtwuxQzLBNrWUnSlN72mT2qqqgAgKNI3qsoKgtA1Meupl2qK28YUYtOD2/36yQHfk8/yTw74AqfNDtIHaAiN9qO2xHvjFDFpwwAoP3KIsLiETrc9GX9PZz0BlOcdQhcYRLBIFIIwYIhE0UtVBc0ABIS2/wilswHf9t5QA0K1mAzee24fnphMYEQk+Vs3MXreOV1u35Ne0yc5bDYKdm0ndfwknyk7LghC18Rvay8117YN3HZUlqOnA74up4xCKbknuF6QJIlRc86iaP8e6stKPHKOQxvWYjG2MmrOWR45viAIniESRS+NPqPtarqpuv03/p4M+MqyTFO1icAw7w7wjpt/PhqdjvXvvIa7+5RYTUa2fPI+scOySBzZ+4q0giD0P/HoqZdiToxNlByqJzo16LSv3zc/65QxCuh4wLehwoix2UbssBCPxdsdfkHBzLjsGta99Qr7Vn3NuHMWuuW4siyz6pUXMDU3sfh3jyBJ3rtz6m9NJhvf59Wyq7iR4gYTRqsDjVJBXIieMQnBzM2KIincz9thCkKnRKLoJf9gLfFZIRzaVEHOOUmofjYLqCcDvrnfl6FQSaSN6/00W3fJWXABRft2s+7NlwgMDyd9wpQ+HU+WZTZ9+DZ5WzYw84rrfpxdNdgV1xv595pjfLGvApvTRYBWRWqEP0F6FRa7iw1Ha/l0dxl/4iDT0sK5+8wMpqWHeztsQWiXaIXaB6VHGvj8X3uZuDCFKYt617uhusjAp0/tYtTseGZf4RtvolaTiU/++ntqigqYd8OtjDlrQa/uAhx2O+vefIn9q79hzJkLOOumOwb93YTLJfPShgL+uSofhQIun5jIReMTGBMfjEJx6msvqTfxZW4Fb28ppspgYfG4OB5dPIogndpL0fs+lyyzt8XED01G8o0WGuwOJCBcoyLTT8eM0ABGB+gH/c9ZR0TPbB+1+s1D5G2rYv6Noxg2oWdTPg31Zpb/fTdIcMUjk9H6+c4bhNVk5Mtnn6Jo7y7SJ05hzrU3EhLT8ZqRnys9lMvaN16irqSIyRdewswrrhsQM51khwPz/v04GxvRDR+OOi6u2/ta7E5+8+FevjlYxbmjYvjzopFEB3U97mSxO/nP98d5ft0xUsL9eOuXk0kIFY+jfsridPFWRR2vltVRammbHRitURGpUSMjU2NzUGtrWyyaptdyU2IkV8eGoRkAP3PuJBKFj7LbnHz+r71UFxmYdlE6485MRFJ0fTVTeayJb185gN3mYvFvc4hMDOyHaHtGdrnYuXIFWz5+D6fdTuaUGQyfNZekkWNQ605/AzQ1N1G4dxe5a7+l/MghAsMjOevG293eBMlTzLkHKP/dvdiLT8z6kiRCLruM6N8/jOJE6ZGOOJwubnlnF2uO1PDIecP51czUHl/Vbi2o5+a3dxKkV7PstulEdSPJDAXbm1q5+0gJRWYbU4P9uTounDlhgURqTr2wqrbaWdNg4L2KenYZTGT4aXl+RDJjA4dO0hWJwofZLA7WvHmYgr21RCYFMuHcZFJGRaBUn3o1I8sytSUt7F9bRt72KoLCdZx76xgiErpfTNAbWhvq2blyBQfWfYfVaERSKAiNjScgLByVWo3NYsZQW4OhtgaA4OgYcuZfwJiz5g+YUh22oiIKL7kUZVAQUb+7F3VCAoaVK2l4622CL1xE3P/9X6f7/+2rw7y8oYDHFo/i2qnJvY5jf1kTV7y8layYQJbeMg21jzW36m8fVTZwb14J8VoNT2clMjus6wsqWZZZVW/gwfwy6mwOnh+RzKKoEM8H6wNEovBxsiyTv72abZ8X0FJvQa1VEpUcSGCYDkkpYW6xU1faQmujFZVaweg5CUw8LwWNbuDMJ3A67JQezKXs8EHqy0owNjXgdDhQa3UEhkcQmZxK0qixRKcNG3DPiEtvux3Tzp2kfbbilMdNtf9+jroXXyTp7bfwnzy53X23FdRz+ctbuXpKEo9fNLrPsazcX8kd7+/mt2dncveZGX0+3kC1vLqR2w4Vc0ZoIK+MSiFI1bP+JQ12BzfkFrKj2cibo1M5J6L9KgqDiUgUA4TL6aLkUAPFB+qpLWnB2GTF5ZLR+asJi/MnISuU9PFR6Px9ZzxiqLNX13BszhwibruNzbOWnDJT7f65qQy/+0r8Z88m/umnTtvX5ZI5/7lNNJvtrP7tGeg17mnGdMf7u1l9qJr1980lJnhg3JW5U57RwvydeYwL9OPDsenoenlnZXQ4WbL3GMdNVtZMyiJZr3VzpL7FU4li4FzODhAKpYKU0RGkjPZMO1HB/cx794IssytuxGn1uR78Mo83R+Sg2LWr3X3X5dVwqNLAPy4b67YkAfDggmy+OVDFqxsLeOT8EW477kAgyzIP5JXip1Tw8siUXicJAH+VkldHpTJ3+xEeOVrOO2N6NztxqBvaD0AFAWha9ikAL+xtaLc+1+YWNY4OGi19sL2EqEAtF4zt/uyo7kgM82PBqBg+2V2GzdF+mZjBalNjK1ubjdyfGkuUtu933ok6Db9OjmZVvYHdBqMbIhx6RKIQhjx1VFvXvabG9vuBOIxGJO3pjyxMNgcb8uu4YGycRwadLxoXT5PJzrbCercf25e9W1lPmFrJFTFhbjvmDfER+CkVvFcxtL6X7iIShTDkBS+5CIAcZ/tvIiNaKtCmp5/2+d3FTdicLs7I9MyK+unDwlEpJH44PnTe3JyyzNp6A+dGBHf5yKm7/V4AAlRKzg4PYnW9we11zIYCkSiEIU8/ejSqyEiuq9mFXnXqr8TI1grS64oIOHPeafsdqmwrNT86vnuzaXryxgbgp1ExLCqAw5WGbr6Sga/QbKXF6WJycOdTxk/2eylvMqPDAk0lPLdsDV/syO9wn0nB/lTbHFTbuu7iKJxKDGYLQ56kUhF+6y04HvsrLw8fw0PBkyhvtjJSbeGvhz5GGR5O6OWXn7ZfRZOFQK2KUP/OF+NBzxpZ/VRimB/F9UPnuXq5xQ5Asr6T76nTwZav3uVxNjBZe4QE6SfjRyuBDXGQPA1GXAhZC0HZNs6RrNOcOIeNGDeMfQwlIlEIAhB65ZVY9u+HT9/mrdSNqOPjMe3ahSRJJLz8Esqg0ysEGyx2gvTde8PpaSOrk4L1alosQ+cKeEtT2zhRXXtX/bIMB5fBmkd5yl5EgyKALa5RvO+aRx3ByEhEYOCBZAcUrocDn0JQAsz7PYy9Ev2JR1lWl3j01FMiUQgCICkUxD75JP4zZtD82ec4GxoIvuACwm/8FZqkpHb3Wba780dHP9V+wyoXVebjfJLfQHlrOUa7EbVCTaQ+ksywTMZFjkOWQTHAFi/2xfSQAJ4trj59fMJigM9uh8NfQMxoHlbfz9KWUTh+9hYWH6LngUvmgcsJR1fBhqdgxW1w4FOa5zwLQKBKPHHvKZEoBOEESZIIXrSI4EWLurX9DTNS+GhHKbIsd7kSPS5ET/mJZCEpW1GHbUEdvBOF2sBffgCVpMJP7YfdZcfsaNtOo9AQ4BqL1v+Mvr2wASTLv21xYZHZ+r9Pmhrg7UVQcxjOfhSm3cnkfVUsX5aLo6N+LwolZC2AjHNg52vw7cPkq16F6MWkDPJFd54gEoUg9FJiqB8mm5PaVitRgZ2vnm5rZLUPR8AmtJHfgcKGbMpmSepN3Dj5TOID4lFIbVe6zdZmDtYdZH3Zej44tAI5aAf3rT/AQ1MeIkznvimjvihGqyZRp2FDQws3JkSCwwrvXw61+XDlR5DR1ka32/1eFAqYfBNEjWDD3uMMt5QT6MoA/Pv5lQ1sIlF4mCzLOGpqcLW2ogwJQRkWNuDqIAntGxnXNm6xr7SZs0d0nijmjw7j/aLPOGzYgqM1k2DTJTx41ux2xyeCtcFMj5/OyLBJvPxZFrMmHmJNyTJ21+zm+XnPMzx8uEdej69YGBnM62V1VFvtRG94HMq2w6Vv/pgkTlqcE9/p+M5PFUVPZGuwH/cWvQnfrYbz/+n+wAcx8bDOQ2xFRVT+4Y8cnTadY2fMoeC88zk6YybH5p1Jzd//jqN+6MyNH6zGJoagVytZl1fT6XY2p4271tzFEcMPPDDpAQ7d/glb77uyyze5jUdrkV1a7sy5nffPex+lpORX3/6KIw1H3PkyfM4N8RG4kPnXkXzY8jzkXAsjL+rTMZ8pqkKjUHBNbDjsfB0q9ron2CFCJAoPaHj/fQouWETzF1/gP3sWMX/6I3F//zvRDz2Ibvhw6t94k+PnLqRl3Tpvhyr0gU6t5MzhUazcX4nlZzOafurJ7U+yrWobj898nGtGXNPtO8pPd5URFahlQnIo2WHZvLXgLfzUfty19i6arc3uehk+J0Wv5dq4CN6qt7I1NAfO/GOfjre23sDHVY3cnBBJzJxfgy4YNok7ip4QicLNGt5/n+pHH8N/+nSGrfqO+KeeIvTKKwk+/zzCfvELEl98gbQvPkeTkEDZnXfRummzt0MW+uCaqck0m+18tKO03a//UPEDH+d/zA0jb+CC9Au6fdy8qhbW5dVy5eQklCcaYcUGxPLs3GepM9Xx9I6n3RK/r/p9YijJlkp+NfKvHFecPjW5uw61mrntUDEj/HX8NiWmLUmMuwaOrARzk/sCHuREonAjW2kpNU88ScAZZ5DwwvOoItsv7aBNSyPp7bfQpqVR+dBDuEymfo5UcJcpqWFMSQ3j32uO0mSynfI1WZb51+5/ER8Qzx05d3T7mLIs89eVhwjUqbhhRsopXxsZMZJrR17L58c/53jTcXe8BJ8UWLmTd3IfAJWWC3cfY1tT+3W4oOMV7+sbWliy5xh+SgVvjk79cR0Fwy8Alx2KNvbHSxkURKJwo8Z33wMg5tFHkVRt8wQ6+iFWBgQQ8+c/4aitpfnLL70Ws9A3kiTxpwtG0my289Cy3FPqCO2v28+h+kP8ctQv0Sq7PyXzna3FbDxax/3zswjxO32F8g0jb0CtUPNR3kdueQ0+qSqXYeZSPhsVT6BKwUV7jvHI0TJqbfZTNvtpKQ+ZthXvD3x5kEs3H+HyfceJ1qpZkTOMpJ9OiY3LAUkBVbn9+5oGMJEo3Mi4dSt+U6agjo4C2v8hfmhZ7o/JQj9+PKq4WExbt3oxaqGvRsQFcf+CLL4+UMUz3/2v1tCakjWoFCrOTT23y2OcvKBIeXAlf/zsICNig7h6SvstVUN1oZyReAZrStYM3gJ3LZWg9mNYeDTfTczimrhwXi+rY+IPh7jpQBHvV9az12DiiXVHMSnBFaTGkeCHbVwYzdMj2Wgxc0tCJF9NyDi9WZFaB/6RYKjwzmsbgMT0WDey5uVhzcv78d9dlW2QJAl1XByO2vZ7HQgDx02z0iioNfL8umOY7U4eOjeb/bX7GRE+gkBN532ef14HCqCgtpXP91V0ODNqUswkVhWvotpUTYx/jFtfi0/44fkfPwxUKXkqK5FbEiN5vayOz2ub+KK2qe2LOSGn7mdxoixpRVVi5C9/yun4+EotuIZOaZS+EonCzX7ab7n9sg2nft5ZV482Y5jH4xI8S5Ik/nbRaHRqJa9tKmRPSSON4WVMj5/S5b5Pfn3ktAsKi8PVaR2o5KC2u42ylrLBmShm3wcbnganA5Rtb1Ppfjoez0zgsYx4jpusHDNZuP/zgzSabEg2F5LRjmRyItFWyqNDsgzmRtD2fpB8qBGPntwo5NJLcDY14WxpAdrKNrTn5OetR49iKypCP2FCv8UoeI5CIfHnRSN59opxFNQZqbNU8/Xx7/k6t/K0gW6L3cnWgnoeWZFLlcHS7vE6utAA8FP5tR3H2f6+A17Yif4f9UdP+5JCksjw13FuZAiPjU8hsMaKstaC4kSSOKWUR3sMFWBrgXBxgdZd4o7CjUKvvJKmTz6l5qmniXn0LyfKNpz6SOHkD7Fst1P12F9R+PsTfOGFXoxacLcLx8VzRmYksz4GhzWc297bDUBEgIYgnRqz3UlNixWnS0arUuCnVmJqZx1GRxcaACZ720w5vaqTK+eBLGlq298F30NUxyvRu13K46cK1p16DqFLIlG4kW7ECMJv/BX1r7yKpNWy6L7fAaNP+yG+INWfsl//BtP27cQ+8QSq0FBvhy64WYifhkkxkzDZTdxz1lR2lTRS2mCixeJAp1YSF6xjVHww04dFsPpQdYcXFB0pNBQCkBCQ4PHX4hVhqRA9Gva+B1NuhU4WKfaklAcAez+AkGSIGe2GQIeGfk0UkiS9BowAVsqy/Nf+PHd/ibznHmSbjYa33qZl7RpmLLmYs2fnoAiKxllfh/G79zi2bBkuo5HoRx4h5KLF3g5Z8JBxkeN448AbZMermJLW8WOO3lwVb6vcRqx/LFF+UW6P22dMvgm+uBvyvobshe45ZuFGKN4E5zzeafIRTtVviUKSpCWAUpblaZIkvS5JUoYsy6c/gBzgJIWC6IceImDOHOr+81/qnn/+1A1UKgLmnEHknXeiy872TpBCvzgz6UxeyX2FlQUruWr4VZ1u25Or4jpzHZvKN3F51uWDu8DkuKvghxfg6/sheTroQ/p2PJsRvrwHgpNg0q/cEuJQ0Z93FHOApSc+/g6YCQy6RHGS/7Rp+E+bhqOxEWv+UVzGtuqx2swslAGixPFQMCJ8BGMix/DGwTe4KOMit40nvJb7Gk7ZyeVZp7dnHVSUalj8Irw+Hz65Aa78EFS97CXhdMCym6HhOFy7HNSDdGzHQ/pz1pM/cLIlWAMQ/fMNJEm6WZKknZIk7aytre3H0DxHFRqK/5TJBM6bh9/48SJJDCGSJHHP+HuoMlbxr13/cssx99bs5f0j77MkYwkpwSluOaZPS5gI5/8Ljq9t60thbuz5MawtsPQ6OPIlzP8bpM1xd5SDXn8milbgZBoPaO/csiy/LMvyRFmWJ0Z2UCdJEAaSiTETuWb4Nbx/5H2W5i3teodOlBj+v727j83qLOM4/v1Bi7zJVt5FbRGELeDAQjXgkGCQLSwkbv4zI910uswsmlk7M12QxGzZH2aiWfbiNDW+sBeDJlaDW8AtM8NFdNSIY5rNJitjiwh1Dhg6YA+Xf5xT9gTa026hz2nP+X2SJs9zn3Oaq3eu9uo55znX/SJtT7QxZ9Ic2pe3n6cIR4Fl18An7oOe38P3ViUN/YbyRHoEdD8G938Unn8U1t8JK24c/ngLqJaXnrpILjftBpYCz2XvblYM7S3t7D+6n9t3386rJ17l+kuuP7Oa3VDtPbyXtifaqESFe9beM+jT3oXTvBFmXJysm/2zT8PsJdDcCvPXwtR5yUp2kBSH/7yQnIH85SF4uSvZ/pntMPfSfH+GUUy16hUjaQqwC3gcWA+siIgBm+q3tLTEnj17ahKb2XA7WTnJ5qc288gLj7B81nJu+dAtLJq2aNDjjpw4QsczHWz921ZmT5rNvWvvZf6F82sQ8QhVOQV7H4Y/fh/+tS8Zq5uQ9G4ScLwX0mdMmH5R8smpZde+/Xsbo4ykrohoOe/ft5ZNxSQ1AOuAJyPiYNa+LhRWNBFBZ3cnW7q2cOTEEVpmtXD53MtpntlM05QmxteNp3K6wqH/HuLZfydrZu/o2cHrb7zOVQuu4uaWm5kyzm0nzujthv1PQe/zSYEgYOI0mL4AGj8CMy4q3UdgC1Eo3goXCiuqYyePse25bXR2d9JztOfMeN2YOiqnKwTJ7+Tk+smsa1pH66JWFjYszClaG01cKMwKJiI4cOwA+3r38dJrL3H81HHqx9Qzc+JMFjYsZPG0xdSPrc87TBtFhqtQuIWHWU4k0TilkcYpjXmHYpbJ3WPNzCyTC4WZmWVyoTAzs0wuFGZmlsmFwsysBCRNlbQqfX2DpK+lr6+QlPn0pwuFmVl5dEh6P0lj1oqkucBdwIBdMsCFwsysFCLiFeAG4IKq4WnAdRHxcv9HJfwchZlZwUn6ArAROJ0OzQDGAxvS7XXAAxFxf3/Hu1CYmRXfj4GOiKhI+gpwBXAh8KuI+K6kd5C0VeyXC4WZWcFFxAkASU3AdcAWkrOKj0n6R0Rszzre9yjMzEpAUjPwG+CrwPF0+LPAZkk/lbR4oGNdKMzMCi4tAr8AvhQRO/vG0xvcq4FukjOMfvnSk5lZwUXEs5IujohT6dA40r//6WWp27KOH7FtxiUdBvbnHUeV6UBv3kGMAJ6HhOch4XlIjJR5aIqIAc8M3q4RWyhGGkl7hqPP+2jjeUh4HhKeh0TR58H3KMzMLJMLhZmZZXKhGLof5B3ACOF5SHgeEp6HRKHnwfcozMwsk88ozMwskwvFEEj6oaQ/SPpG3rHkRVKdpBcl/S79uiTvmGpN0ixJu6relzIvquehjHkh6QJJj0raKemXksYVPRdcKAYh6ZPA2IhYCcyTtCDvmHKyBHg4ItakX8/kHVAtSWoAfgJMSt+XMi/OngfKmRcbge9ExGXAQeBTFDwXXCgGtwbYlr7eCazKL5RcrQA2SPpT+t9T2Z7qrwBXA0fT92soZ16cPQ+ly4uIuC8ifpu+nQG0UvBccKEY3CSgb1GPV4BZOcaSp6eBj0fEh4F6kjbFpRERRyOiehWwUuZFP/NQ2ryQtBJoAA5Q8FxwoRjca8CE9PVkyjtnf42If6av9wCFO71+i5wXiVLmhaSpwN3A5yhBLhTuBxoGXbx5KrkU6MkvlFxtlbRU0ljgSmBvzvHkzXmRKF1eSBoH/By4NSL2U4JcKPz1xPOgE9glaQ6wnuSabBndBjxEsgrWryPisZzjyVsnzgsoZ158HlgGbJK0CfgRcE2Rc8EP3A1B+kmPdcCTEXEw73hsZHBeWJ+i54ILhZmZZfI9CjMzy+RCYWZmmVwozIYobdXw7gG2zat1PGa14kJhNnR3Au8dYNuNklbXMhizWnGhMKsi6ZuS/l7V5O6D6fhK4FBE7B7g0FuBmyT5d8oKx89RmJ3rjoh44Kyxa4GvD3RARLwh6XHgUmDXQPuZjUb+78dsaCb09TiStFHStyVNlvRnSe9M99kNNOcXotnwcKEwO9emqktPY8/eGBEPAh8A7gK+FRHH0k3/482eP2aF4UJhdq47qtZXqKRjYySpap8OYANJz58+7yPpJGpWKC4UZkPzNLAaklXdgJuArcAXq/a5kmQ9ArNCcaEwO1f1paer07EO4MuSJgJtwHaSTzq1SnqPpFXA4YjozSdks+HjXk9mQyRpCfCuiNjRz7Z24O6IOFX7yMyGlwuFmZll8qUnMzPL5EJhZmaZXCjMzCyTC4WZmWVyoTAzs0wuFGZmlun/tIaRz9muOv4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "if __name__== '__main__':\n",
    "    file_path = './使用数据'\n",
    "    out_path = './输出数据'\n",
    "    csv_file = 'sample_plot.txt'\n",
    "    infile = os.path.join(file_path,csv_file)\n",
    "    plotID = csv_file.split('.')[0]  # get from csv\n",
    "    outfile = os.path.join(out_path,plotID+'_2.shp')\n",
    "    canopy_crow_create(infile,outfile,20)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
